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Assembly  and  Application  of  a  Shock  Data  Analysis  System 


1.  Introduction 

Measured  data  from  shipboard  shock  environments  usually  consist  of 
analogue  magnetic*tape  recordings  of  acceleration  and  strain  versus  time. 
Usually,  little  analysis  is  required  of  strain  recordings  save  spectral  analy¬ 
sis  and  perhaps  sums  and  differences.  From  the  accelerometer  recordings, 
however,  it  is  usually  necessary  to  extract  velocities  and  displacements,  and 
often  shock  spectra,  as  well.  Shock  data  analysis  facilities  tend  to  run  to 
two  extremes  of  size  and  capacity.  At  one  end  is  the  large  digital  processing 
facility  where  hundreds  of  records  (from  an  at-sea  ship  shock  test,  for 
example)  are  processed  in  batch,  usually  after  the  test  is  completed.  At  the 
under  end  is  the  small,  often  analogue,  facility  which  can  do  limited  analysis 
of  a  small  number  of  records,  although  perhaps  very  slowly.  The  latter  may  be 
used  to  advantage  for  analyzing  a  few  selected  critical  records  in  the  field, 
or  in  the  shock  simulation  lab. 

A  need  exists  for  a  capability  in  the  middle  ground  -  a  system  capable  of 
doing  analysis  on  a  modest  number  of  records  in  a  reasonable  time,  which 
accepts  analogue  magnetic  tapes  as  input  and  produces  report-ready  graphs  and 
tables  as  output,  and  which  is  field-transportable  if  necessary. 

This  report  describes  such  a  system  which  has  been  assembled  with 
standard,  off-the-shelf  instruments  which  are  marketed  nationwide  and  easily 
available.  It  is  not  in  any  way  optimized  for  the  specific  task  of  shock  data 
analysis,  hence  may  not  be  competitive  in  speed,  cost,  etc.  to  a  system  which 
has  been  optimized  in  hardware  and/or  software  to  this  specific  application. 

It  furnishes  an  example  of  what  can  be  done  along  these  lines  with  virtually 
no  expenditure  of  R  &  D  effort,  and  is  an  example  in  that  any  of  its  component 
instruments  could  be  replaced  by  a  different  supplier's  product  with  little 
change  in  operation  or  performance. 

2.  Background 

2.1.  Shock  motions  may  be  described  in  terms  of  acceleration  -, 
velocity  -,  or  displacement-time  histories,  and  in  all  three  parameters  may 
vary  greatly  in  waveform  complexity.  Simpler  shock  motions  may  be  closely 
defined  by  a  few  waveform  parameters,  and  these  may  also  provide  a  reliable 
basis  for  comparison  of  shock  motions.  Shipboard  shock,  however,  is  extremely 
complex.  Motions  are  typified  by  waveforms  of  high  frequency  range,  high 
dynamic  range  and  long  duration  which  cannot  be  defined  adequately  by  a  few 
parameters.  The  additional  descriptions  found  most  useful  for  a  motion  of 
this  type  are  the  shock  spectrum,  which  may  be  interpreted  as  a  measure  of  the 
effects  of  a  shock  motion,  and  the  design  shock  spectrum,  which  provides  the 
data  needed  to  design  a  structure  to  withstand  it.  Each  provides  a  basis  for 
comparing  shock  severities,  and  applies  equally  well  to  shock  motions  of  all 
degrees  of  complexity  of  waveform.  The  general  nature  of  shock  motions,  shock 
spectra  and  design  shock  spectra  are  discussed  in  Refs.  (1)  and  (2). 

2.2.  While  a  shock  motion  may  be  identified  equally  well  by  its  acceler- 
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ation,  velocity  or  displacement  waveform,  one  of  them  will  usually  be  most 
convenient  for  a  given  application.  However,  in  most  situations  only  acceler¬ 
ation  can  be  measured  properly.  This  is  almost  always  the  situation  in  the 
case  of  shipboard  shock,  and  the  bulk  of  data  in  this  field  will  consist  of 
records  of  acceleration  and  strain  versus  time,  with  perhaps  a  small  fraction 
being  records  of  velocity  versus  time.  Thus,  an  analysis  system  for  shipboard 
shock  motion  records  must  be  able  to  produce  time  histories  of  acceleration, 
velocity  and  displacement,  as  well  as  shock  spectra  and  design  shock  spectra, 
given  input  signals  representing  acceleration  -  or  velocity-time  histories, 
and  preferably  from  displacement-time  histories  as  well.  Since  shipboard 
shock  motions  are  among  the  most  demanding  in  terms  of  functional 
requirements,  an  analysis  system  capable  of  handling  them  will  also  be  able  to 
handle  shock  motions  of  most  other  types. 

3.  Analysis  System  Characteristics 

3.1.  The  principal  application  of  the  system  is  the  analysis  of  shock 
motions  measured  on  laboratory  shock  machines.  This  includes  those  which 
generate  simple-pulse  shock,  but  most  prominent  are  those  which  generate 
shipboard  shock  type  environments:  the  Lightweight  Shock  Machine  (LWSM),  the 
Mediumweight  Shock  Machine  (MWSM) ,  the  Floating  Shock  Platform  (FSP)  and  the 
Large  Floating  Shock  Platform  (LFSP)  (Ref.  3).  Tests  on  these  machines  will 
be  lightly  instrumented  compared  to  at-sea  tests  of  ships  and  submarines,  or 
some  of  the  more  specialized  test  vehicles.  The  typical  instrumentation  suite 
may  be  a  few  channels  on  the  simpler-pulse  machines  and  the  LWSM,  up  to  a 
dozen  on  the  MWSM  and  FSP  and  a  few  dozen  on  the  LFSP,  while  on  some  occasions 
several  times  as  many  channels  may  be  required.  Normally,  about  one- third  of 
the  total  channels  will  fall  in  the  category  of  equipment  input  motions,  while 
most  of  the  remainder  will  be  equipment  response  motions  and  structural 
strains. 

3.2.  The  function  of  the  analysis  system  is  to  accept  analogue  input 
signals,  perform  whatever  processing  is  required,  and  produce  output  represen¬ 
tations  of  the  desired  quantities,  scaled  in  the  appropriate  engineering 
units,  in  graphical  or  tabular  form,  and  in  video  or  hard-copy  display.  All 
signals  require  A-D  conversion  and  storage,  scaling  and  display.  All  motion 
signals  require  in  addition  computation  and  display  of  acceleration,  velocity 
and  displacement,  each  properly  scaled,  which  requires  the  functions  of 
integration  and  possibly  differentiation.  Input  motion  signals,  and  some 
response  motion  signals,  also  require  the  computation  of  maximax  and  residual 
shock  spectra.  This  last  computation  is  by  far  the  most  complex  and  time- 
consuming,  but  is  needed  for  only  a  few  channels.  For  a  suite  of  a  dozen 
channels  of  assorted  information,  the  system  will  perform  a  complete  analysis 
in  an  hour  or  two. 


Shock  spectrum  analyzers  are  available  commercially,  and  are  typically 
compact  units  and  very  fast.  As  stand-alone  systems,  they  are  lacking  in  that 
they  cannot  do  a  complete  analysis  (integration  and  scaling,  etc.)  and  have  no 
display.  As  components  for  calculating  shock  spectra  they  are  wanting  in  that 
they  cannot  compute  undamped  spectra  and  their  analysis  frequencies  are  fixed. 


3.3#  Since  Che  analysis  system  is  an  assembly  of  component  units  designed 
for  use  as  independent,  general  purpose  instruments,  it  has  a  great  degree  of 
flexibility.  With  different  switch  settings,  it  would  function  equally  well 
for  analyzing  earthquake  motions  or  blast  loadings.  A  brief  description  of 
each  functional  block  of  the  system  and  its  capabilities  follows. 

3.3.1.  Magnetic-Tape  Recorder 

The  primitive  form  of  shock  data  records  is  usually  FM 
magnetic  tape.  The  recorder  used  with  the  analysis  system  has  14- 
track  record/reproduce  capability,  plus  edge  tracks  for  voice  annota¬ 
tion.  The  normal  configuration  has  two  housekeeping  channels,  one  of 
them  DR  for  a  tape-speed  servo  reference  signal  and  Che  other  an  FM 
channel  devoted  to  time-code  in  IRIG  B  format,  and  12  FM  data 
channels.  All  FM  channels  are  IRIG  Intermediate  Band  I.  Normal 
recording  speed  is  30  ips,  giving  a  nominal  bandwidth  of  0  -  10 
KHz.  The  recorder  can  be  operated  at  tapes  speeds  from  15/16  to  120 
ips,  allowing  the  traditional  analysis  tricks  of  data  speed-up  and 
slow-down;  record  density  can  be  changed  by  rearranging  plug-in 
components  on  the  electronics  cards.  Output  level  is  normally  set 
for  rras  for  i40%  peak-to-peak  (sinusoidal)  modulation.  This 
particular  unit  was  selected  because  it  performs  a  set  of  canned 
calibration  and  set-up  routines  under  control  of  a  built-in  micropro¬ 
cessor,  and  because  it  is  a  portable  unit.  Both  features  are 
convenient  for  field  operations. 

3.3.2.  Auxiliary  Equipment 

The  principal  auxiliary  equipment  units  are  a  time-code  gener¬ 
ator/reader,  a  direct-writing  oscillograph,  and  a  storage  CRO.  In 
generator  mode,  the  time-code  unit  produces  an  IRIG  B  code  for 
recording  in  the  magnetic  tape.  In  reader  mode,  it  monitors  the 
recorded  time-code  and  produces  a  pulse  when  a  chosen  time  value 
occurs.  This  time  can  be  selected  to  the  millisecond  and  provides  a 
start  signal  to  the  transient-storage  unit  (discussed  below)  and 
facilitates  synchronization  of  the  recorded  data  signals.  The  other 
units,  the  CRO  and  oscillograph,  provide  quick-look  display  of  the 
recorded  data  signals  and  time-code.  The  CRO  is  a  two-channel, 
single  time-base  type  similar  to  many  currently  available.  The 
oscillograph  produces  strip  charts  of  up  to  12  channels  with  a  band¬ 
width  of  0  -  5  KHz  at  paper  speeds  from  0.1  to  160  ips. 

3.3.3.  Data  Filter 

The  data  signals  output  from  the  tape  recorder  are  passed 
through  an  8-channel  data  filter.  Each  channel  has  a  low-pass 
characteristic  with  cutoff  frequency  selectable  from  1  Hz  to  15  KHz, 
and  is  direct-coupled.  The  unique  features  of  this  unit  are  that  the 
filter  characteristic  is  optimized  for  linear  phase- shift  within 
pass -band,  minimizing  envelope  distortion,  and  that  phase-shift 
between  channels  is  closely  controlled  at  no  more  than  2®  between  any 
pair.  The  filter  characteristic  also  features  a  cutoff  rate  of  43 
db/ octave ,  so  that  aliasing  is  negligible  for  sampling  frequencies  of 
4,2  times  the  cutoff  frequency  or  higher.  Each  channel  has  unity 
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gain,  and  can  accept  inputs  of  up  to  10  «  As  used  in  the  analysis 

system,  the  cutoff  frequency  is  normally  set  to  IKHz. 

3*3*4.  Transient  Storage  Unit 

The  transient-storage  unit  contains  seven  independent  AD  and 
digital  storage  channels  with  a  common  time-base.  Each  channel  can 
be  operated  in  an  inhibit  mode  which  preserves  the  existing  stored 
signal.  The  unit  can  be  operated  in  a  self-triggered  mode,  where  the 
signal  in  any  channel  can  be  used  to  initiate  the  capture  of  the 
signals  in  all  channels  which  are  not  inhibited.  As  used  in  the 
analysis  system,  the  capture-initiate  signal  is  derived  from  the 
time-code  reader  to  preserve  commonality  of  time-base  between  succes¬ 
sive  batches  of  signals.  Upon  initiation,  each  channel  captures  4096 
10-bit  samples  of  its  signal,  for  a  dynamic  range  of  60  db.  As  used 
in  the  analysis  system,  the  channel  full-scale  range  is  ±  2  ^  ;  it 
can  be  varied  from  ±0.1  to  irS^.  The  sampling  period  used  is  100 
microsec.,  for  a  total  record  length  of  400  mlllisec.;  the  sampling 
period  can  be  varied  from  5  mlcrosec.  to  2.5  sec.  Each  channel 
provides  outputs  of  the  digitized  "staircase"  and  smoothed  analogue 
versions  of  the  stored  signal  waveform.  As  used  in  the  analysis 
system,  the  signal  stored  in  any  selected  channel  is  read  out  on  the 
IEEE-488  bus  in  bit-parallel,  character-serial  ASCII  code.  This 
particular  unit  was  selected  because  of  its  7-channel  capacity  and 
the  fact  that  each  A-D  converter  and  storage  channel  is  functionally 
independent  from  the  others. 

3.3.5.  Processor  and  Display  Unit 

The  processor  and  display  unit  is  a  large-screen  desk-top 
computer  and  peripherals.  The  computer  has  64K-byte  memory  and  a 
resident  BASIC  interpreter,  leaving  some  53K-bytes  of  usable 
memory.  Additional  ROM  packs  provide  signal  processing  routines, 
including  integration  and  differentiation,  and  text-editing 
capability.  The  principal  display  is  a  built-in  large  screen  CRT. 

The  peripherals  include  a  hard-copy  unit  which  gives  a  permanent 
record  of  the  CRT  display  with  reduced  size,  and  a  digital  X-Y 
plotter  which  gives  high  quality  graphs  up  to  11  x  17  in.  Other 
peripherals  are  a  high  capacity  disc-file  manager,  which  allows  a 
large  library  of  programs  to  be  resident  on  floppy  disc,  and  an 
additional  tape-cassette  drive-unit,  allowing  more  flexibility  in 
data  transfer.  All  of  these  components  are  produced  by  the  same 
manufacturer,  and  are  configured  and  interfaced  to  form  a  compatible 
system.  This  (processor  and  display)  unit  was  selected  because  it 
possesses  the  requisite  computational  power  and  has  an  excellent 
repertoire  of  graphics  routines.  Other  features  are  good  man-machine 
interface,  BASIC  interpreter  and  large-screen  display.  These  seem  to 
assist  the  Intended  user  profile,  engineering/technical  personnel 
without  extensive  familarity  with  computing  systems,  in  learning  to 
operate  this  analysis  system  with  a  minimxm  of  instruction.  The 
existence  of  a  comprehensive  assortment  of  proven  peripherals  was 
also  a  factor  of  consideration.  Another  feature  of  the  computer  is 
that  it  can  also  be  operated  as  a  high-speed  (9600  baud)  intelligent 


4 


terminal,  should  operation  with  a  larger  host  computer  be 
desirable*  The  last  peripheral  in  the  system  is  a  high-speed 
printing  terminal  from  another  manufacturer.  When  used  with  the 
analyeiit. system,  it  functions  as  a  moderately  fast  line-printer* 

3*4*  Ttie  functional  block  diagram  of  the  analysis  system  is  in  shown  in 
Fig*  1* 

4.  Computation  of  Shock  Spectra 

4.1*  Recalling  that  the  shock  spectrum  of  a  shock  motion  may  be 
interpreted  as  the  maximum  relative  displacements  of  a  set  of  massless  single- 
degree-of-freedom  systems  whose  bases  are  forced  to  follow  the  notion,  it  is 
apparent  that  shock  spectra  may  be  derived  for  any  degree  of  damping*  For  the 
purpose  of  comparing  shock  motion  environments  and  specifying  design  shock 
spectra,  the  undamped  version  is  preferred,  since  the  damped  spectra  are  less 
sensitive  to  variations  in  the  input  motion.  Shock  spectra  can  be  computed 
from  the  differential  equation  of  motion,  or  from  the  corresponding  integral 
equation*  Either  form  can  be  expressed  to  involve  the  input  motion  in  terms 
of  acceleration,  displacement  or  velocity.  The  procedure  used  here  is  that  of 
numerical  solution  of  the  differential  equation  relating  relative  displacement 
to  input  accleration*  The  computational  procedure  requires  the  input  to  the 
program  as  the  acceleration-time  history  of  the  shock  motion.  The 
computational  procedure  uses  the  algorithm  described  in  Ref.  (4). 

The  equation  of  motion  for  an  undamped  oscillator  of  frequency  o  is 
( 1)  X  +  (ii^X  =  -Z 


with  the  solution  for  the  relative  displacement  X 
(2)  X(t)  =  Xq  cos  wt  + sin  wt  -  /qZ(t)  sin  u)(t  -  t)  dt . 

The  numerical  integration  equation  derived  by  G*  O'Hara  (Ref.  4)  is 


(3) 

^n+1^  =  X^u)  cos  a)h  +  X^  sin  wh 


^n(  1  -  cos  mh)  _  cz  /I  +  cos  a)h 
wh  ^n- 1  2u)h 


sin  u)h\ 


with  a  parabolic  approximation  for  the  digitized  input  function 


Z  (t); 


and 

(note  that 
so  chat 


Sn  -  Sn_i  -  Zn+1  "  2Zn  +  Vl 

is  the  second  difference,  not  a  square) 
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In  these  equations,  h  is  the  time  interval  of  integration. 

By  eliminating  in  (3),  R.  Bort's  implementation  of  equation  3  (Ref.  5) 
is 

(5)  X(t  +  h)  =  2  X(t)  cos  wh  -  X(t-h)  +  h  [Z  (t  +  h/2)  -  Z(t  -  h/2)] 

Equation  (5)  has  been  applied  in  this  program  because  it  consists  of  only 
three  additions.  As  implemented  here,  the  term  ^n  brackets 
[Z  (t  +  h/2)  -  Z  (t  -  h/2)]  is  identified  as  hZ  (t),  a  sample  point  from  the 
digitized  input  acceleration  record.  Before  processing  input  data  the 
2  cos  u)h  term  is  determined  and  stored  in  array  C.  All  scaling  factors  are 
deferred  to  the  end  of  the  calculation  for  each  analysis  frequency. 

In  order  to  cut  down  calculation  time,  the  fact  that  lower  analysis 
frequencies  require  fewer  periods  to  reach  peak  response  than  higher 
frequencies  is  employed.  The  number  of  periods  per  analysis  frequency 
considered  in  the  processing  routine  is  given  by 

^  °  Yllg  T 

where  F  »  a)/2iT  . 

The  second  procedure  to  speed  up  the  program  selects  data  points  with  respect 
to  cutoff  frequency  and  analysis  frequency.  The  calculation  time  has  been 
reduced  this  way  from  about  50  min.  to  8  min.  while  producing  only  a  minor 
error  (see  Fig.  2A,B). 

5.  System  Operator's  Manual 

5.1.1.  It  is  assumed  that  a  time  record  of  the  shock  motion  - 
including  a  calibration  step  -  is  available  on  an  analog  tape  and  a 
time  code  is  recorded  on  a  separate  track  of  the  tape.  Playing  back 
the  record  (or  calibration  step)  and  the  time  code  on  the  Visicorder, 
determine  the  millisecond  preceding  the  starting  time  of  the 
record.  Set  the  time  code  generator  to  open  the  time  gate  and 
trigger  the  transient  recorder.  The  transient  recorder  accomplishes 
the  A/D  -  conversion  and  stores  the  digitized  time  record  in  its 
memory .  An  interface  bus  allows  the  computer  to  have  access  to  the 
seven  channels  (memories)  of  the  transient  recorder.  Further  data 
processing  is  done  under  the  control  of  the  program  listed  in  the 
Appendix. 

5.1.2.  This  program  is  written  in  BASIC  and  runs  on  the  desktop 
computer  Tektronix  4054.  The  three  parts  of  the  program  are  stored 
on  the  disc  in  three  files  (a  back  up  tape  cartridge  contains  the 
program  in  binary  code).  Part  I  reads  the  calibration  step  and  the 
time  record  from  the  memory  of  the  transient  recorder,  and  calculates 
shock  spectra.  Part  II  contains  all  output  procedures.  Part  III  is 
basically  the  same  as  Part  I  except  that  it  reads  a  calibrated 
digital  time  record  from  the  external  tape  unit  rather  than  from  the 
transient  recorder.  Both  Part  I  and  Part  III  use  the  output 
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•  CODE  9*42 
29-n*|l-SI  12.49 


(a)  Maximax  sliock  spectrum 


procedures  in  Part  II, 

5.2  Starting  Up 

The  computer  and  at  least  half  of  the  units  connected  with  the 
computer  by  interface  bus  must  be  turned  on.  To  run  the  program,  set  the 
clock  on  the  disc  drive,  mount  the  unit,  load  the  program  and  enter  **RUN’*. 
Example:  Enter:  CALL  "SETTIM**,  ”20  -  MAR  -81  □  8:30” 

Press :  RETURN 

Enter:  CALL  ”M0UNT”,e)  ,  A$ 

Press:  RETURN 

Enter:  OLD  ”  @  SPRGM/PART  I” 

Press:  RETURN 
Enter:  RUN 
Press:  RETURN 

Message  on  the  Screen: 

SPRGM/PART  I  is  loaded  (default  conditions) ! 

The  corresponding  procedure  for  using  binary  tape  is  either: 

Enter:  FIND  2 
Press;  RETURN 
Enter:  CALL  "BOLD" 

Press:  RETURN 
Enter:  RUN 
Press:  RETURN 
or  simply  press:  AUTOLOAD 

While  running  the  program,  the  computer  indicates  a  need  of  communication 
with  the  user  by  a  "beep"  tone  and  a  message  on  the  screen. 

The  symbol  LJ  indicates  one  space. 

5 .3  Default  Conditions 

Default  conditions  are  listed  in  line  290  to  500  in  the  Appendix.  All 
default  parameters  may  be  changed  at  any  time  by  entering  the  desired  values 
from  the  keyboard  (e.g.  entering  F2  *  500  changes  the  upper  analysis  frequency 
to  500  Hz).  The  current  value  of  every  parameter  can  be  checked  by  entering 
the  parameter  and  pressing  RETURN:  the  current  value  of  the  parameter  appears 
on  the  screen. 

The  most  important  parameters  defined  by  default  conditions  are  explained 
below : 

FI,  F2 :  Lower/upper  analysis  frequency:  can  be  any 

positive  value  between  1  and  1000. 

C0:  Number  of  analysis  frequencies  between  FI  and 

F2:  depends  on  memory  available  (Part  I  up  to 
about  I 10,  Part  III  about  160).  C0  *  1  is 
possible  also. 

F8 :  Cut  off  frequency  of  low  pass  filter  (default 

1000  Hz). 

Nl,  N2:  Lower/upper  value  of  y-axis  in  log  display 
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N6 ,  N7 :  Lower/upper  value  of  x-axis  in  log  display 

D:  Primary  address  for  output  devices: 

32  screen 

3  ext.  tape 

4  plotter 
40  printer 

P:  selects  grid  in  log  display: 

0<P<10  sets  grid  line  at  P 
P>10  no  grid 
(P  must  not  be  zero!) 

0:  indicates  the  part  of  the  program  currently 

loaded  in  the  memory 

1.  Part  I 

2.  Part  II 

3.  Part  III 

(Do  not  change  parameter  0  from  the  keyboard) • 
M9 :  Calibration  matrix. 

A2 :  Offset-correction  matrix. 

5.4.  User  Defined  Keys 


The  first  row  of  user  defined  keys  (No.  11  through  15)  are  used  with 
program  Part  I  or  III.  All  other  user  defined  keys  (No.  1  through  10  and  16 
through  20)  belong  to  program  Part  II.  User  Defined  Key  No.  4  and  ^  allow 
switching  from  program  Part  I  or  III  to  Part  II  or  vice  versa.  Pressing  a 
user  defined  key  which  belongs  to  a  program  part  not  currently  loaded  in  the 
memory  causes  an  error  message  to  appear  on  the  screen. 

The  following  describes  each  user  defined  key  and  the  system's 

reaction. 


5.4.1.  First  Row  (No.  11  through  15) 

User  Defined  Key  No  11  triggers  the  calibration  routine  in  program 
Part  I.  After  the  time-record  has  been  defined  as  a  displacement, 
velocity  or  acceleration  record  there  are  two  ways  to  fill  the 
calibration  matrix  M9  with  the  calibration  factors  for  each  channel. 

First  the  calibration  factor  may  be  entered  from  the  keyboard  in  input 
units  per  1  Volt,  or  secondly,  the  calibration  step  may  be  read  from  each 
channel  of  the  transient  recorder.  In  the  latter  case  the  calibration 
level  has  to  be  entered  in  appropriate  units  before  the  computer  reads 
the  content  of  each  channel  and  displays  it  on  the  screen.  The  user  may 
define  the  calibration  level  by  means  of  setting  the  crosshair  first  to 
the  upper  and  then  to  the  lower  calibration  level.  The  calibration 
factor  (i.e.  input  unit/digit)  for  this  particular  channel  will  be 
assigned  to  an  element  of  calibration  matrix  M9.  The  program  starts  this 
procedure  all  over  again  until  the  calibration  steps  of  all  7  channels 
have  been  read.  If  the  program  encounters  a  channel  set  to  ground  or 
calibration  mode  an  appropriate  message  appears  on  the  screen  and  no 
calibration  factor  will  be  defined  for  this  particular  channel.  The 
program  proceeds  with  the  next  channel. 

User  Defined  Key  No.  12  reads  the  time  record  from  a  defined 
channel  of  the  transient  recorder  and  displays  a  plot  of  acceleration. 
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velocity  and  dlsplacement^vs  time  on  the  screen*  Also  each  maximum  and 
minimum  value  is  printed.  When  running  program  Part  I,  the  calibration 
routine  (key  11)  must  be  completed  before  using  Key  Mo*  12*  With  program 
part  III,  User  Defined  Key  No*  12  causes  the  computer  to  read  the  time 
record  (including  other  important  information  like  calibration  factor, 
time  step,  etc*)  from  the  external  tape. 

The  read  routine  has  to  be  completed  for  each  record  before 
calculating  shock  data  by  pushing  User  Defined  Key  No*  13*  The  computer 
then  calculates  shock  data  of  the  time  record  that  has  been  input  by 
means  of  Key  12* 

Using  the  mathematical  algorithm  described  above,  the  program 
computes : 

1*  Analysis  Frequency 

2*  Pseudo  velocity  Maximax-Shock  Spectrum 

3.  Rel*  Velocity  Residual- Shock  Spectrum  (or 
Fourier  Transform) 

4,  Rel*  Displacements 
5  *  Abs .  Acceleration 

6.  Time  when  maximtm  occurs 

7.  Phase 

While  calculating,  the  computer  displays  on  the  screen  the  serial 
number  of  analysis  frequency  (from  1  at  the  lowest  to  C(&  at  the 
highest) ,  the  analysis  frequency,  the  pseudo-velocity,  and  the 
number  of  data  points  used  for  the  calculation  of  shock  spectra  for 
each  analysis  frequency.  The  computing  time  is  also  given  at  the 
end  of  the  procedure.  ^ 

The  number  of  analysis  frequencies  (1  to  160  ),  as  well  as  any 
interesting  range  of  analysis  frequencies  between  1  and  1000  Hz,  can 
be  chosen.  (Note  that  F2  has  to  be  larger  than  F^)* 

The  time  record  and  all  Important  information,  such  as 
calibration  factor,  time  step,  etc.  can  be  stored  on  a  defined  file 
of  the  external  tape  by  pushing  User  Defined  Key  No*  14*  Additional 
analysis  can  be  done  on  this  time  record  by  means  of  part  III  of  the 
program  (refer  to  User  Defined  Key  No.  12).  An  offset  correction  of 
the  time  record  can  be  made  with  the  procedure  of  User  Defined  Key 
No*  15*  The  first  integral  of  the  time  record  is  displayed  on  the 
screen  and  the  user  may  define  the  zero  line  by  means  of  the  cross 
hair.  After  pushing  any  key  the  computer  calculates  an  offset 
correction  for  this  particular  channel  and  considers  that  at  all 
later  operations*  This  procedure  may  be  repeated  several  times. 

*  The  title  or  description  of  the  time  record  is  stored  in  string  W$  and 
appears  on  the  top  of  the  display.  It  may  be  changed  anytime  by  redefining  W  $ 

^Depends  on  memory  available  for  the  calculation  (program  part  III  allows  up 
to  160  analysis  frequencies) . 


12 


However »  the  user  has  to  realize  that  this  procedure  requires  some  actual 
value  to  be  known. 


5.4.2.  Second  Row  (No.  1  through  5) 

The  results  of  the  shock  spectra  computation  (Key  No.  13)  can 
be  saved  by  pushing  User  Defined  Key  No.  1.  This  causes  the 
computer  to  write  all  data  on  a  defined  file  of  the  external  tape. 
These  shock  data  can  be  read  back  with  User  Defined  Key  No.  2. 

A  numerical  output  of  all  analysis  data  can  be  retrieved  by 
pushing  User  Defined  Key  No.  3.  As  for  all  other  output  procedures 
the  primary  address  of  the  selected  output  device  (screen,  printer, 
or  plotter)  has  to  be  specified  here. 

User  Defined  Key  No.  4  makes  the  computer  load  part  II  of  the 
program  while  keeping  all  values  of  all  parameters  and  variables  in 
its  memory. 

5.4.3.  Thl rd  Row  (No .  16  through  20) 

The  third  row  of  user  defined  keys  produces  a  linear  plot  of 
the  shock  spectra.  The  amplitude  scale  is  automatically  adjusted  to 
the  output.  Any  range  of  analysis  frequencies  can  be  selected. 

User  Defined  Key  No.  16  displays  a  plot  of  pseudo-velocity 
vs .  analysis  frequency  (maximax) 

User  Defined  Key  No.  17  displays  a  plot  of  rel.  velocity 
vs.  analysis  frequency  ( Residual) 

User  Defined  Key  No.  18  displays  a  plot  of  acceleration  vs. 
analysis  frequency . 

User  Defined  Key  No.  19  displays  a  plot  of  rel. 
displacements  vs.  analysis  frequency. 

User  Defined  Key  No.  20  displays  a  plot  of  phase  (0  to  180 
degree)  vs.  analysis  frequency. 

5.4.4.  Fourth  Row  (No.  6  through  10) 

A  logarithmic  plot  of  analysis  data  can  be  retrieved  by 
pushing  any  of  the  user  defined  keys  in  the  fourth  row.  User 
Defined  Key  No.  6  displays  the  Maximax  Shock  Spectrum,  No.  7  the 
Residual  Shock  Spectrum,  and  No .  8  both  of  them.  User  Defined  Key 
No.  9  allows  some  analysis  work  on  the  screen  in  determining  modal 
frequencies  and  appropriate  shock  spectrum  values  (the  Design  Shock 
Spectrum)  by  placing  the  two  crosshairs  at  points  of  interest.  The 
thumbwheels  on  the  keyboard  control  the  movement  of  the  crosshairs 
When  a  keyboard  key  is  pressed,  the  location  of  the  crosshair  is 
marked  and  the  appropriate  analysis  frequency  and  shock  spectrum 
value  is  printed  beside  the  graph-plot.  When  the  key  "0**  is 


4r 


except  phase  is  always  displayed  between  0  and  180  degrees. 


pressed,  the  crosshair  disappears  from  the  screen.  (Note:  this 
procedure  works  with  the  screen  only!) 


User  Defined  Key  No.  10  gives  a  plot  of  a  number  of  different 
shock  spectra  stored  on  the  external  tape.  Each  spectrun  is  drawn 
with  a  different  dark  line  and  the  description  of  the  spectrum  is 
printed  wherever  the  crosshairs  have  been  placed. 

6.  Examples  of  Analysis  System  Output. 

6.1.  Figures  3  and  4  show  the  various  output  formats  which  can  be 
selected.  These  were  made  using  the  hard  copy  unit,  which  produces  a 
reduced-size  version  of  the  CRT  display.  Higher  quality  graphs  can  be 
obtained  from  the  plotter. 

6.2.  For  Fig.  3,  the  input  record  was  a  half-sine  acceleration  pulse. 

The  original,  integral  and  double  integral  are  shown  in  Fig.  3A  as  scaled 
acceleration,  velocity  and  displacement.  The  raaximax  shock  spectrum  is 
shown  in  linear  plot,  without  grid-lines,  in  Fig.  3B,  and  in  log-log 
plot,  with  grid-lines,  in  Fig.  3C*  The  maximax  and  residual  are  shown  in 
log-log  plot,  with  factor-of-ten  grid-lines,  in  Fig.  3D.  Finally,  a 
tabular  listing  of  all  computed  data  is  shown  in  Fig.  3E. 

6.3.  The  input  record  for  Fig.  4  was  an  acceleration-time  history 
recorded  on  a  resilient  structure  aboard  the  Floating  Shock  Platform  with 
shot  standoff  of  20  ft.  The  acceleration,  velocity  and  displacement 
waveforms  are  shown  in  Fig.  4A,  utilizing  the  offset  correction  routine 
described  above:  i.e.,  requiring  the  average  velocity  near  the  400 
millisec.  mark  to  be  zero.  The  raaximax  and  residual  spectra  are  given  in 
Fig.  4B,  and  again  in  Fig.  4C,  which  demonstrates  the  extraction  of 
design  shock  spectral  parameters.  Modal  frequencies  and  shock  spectrum 
values  are  lined  up  with  the  CRT  crosshairs,  and  their  values  tabulated. 
The  values  read  are  also  identified  on  the  graph.  The  tabular  listing  of 
computed  data  is  shown  in  Fig.  4D. 
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HALFSINE  100g- 


1 1 'VI#  record  (a)  Acceleration,  velocity  and  displacement  time  histories 

Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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Maximax  shock  spectrum,  linear  plot,  without  grid 
~  Output  for  analysis  of  a  half-sine  acceleration  pulse 


Maicimax 


Maximax  and  residual  shock  spectra^  log'plot,  with  factorof-te 
Fig.  3  ^  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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(e)  Tabular  listing  of  all  computed  data 
Fig.  3  ^  Output  for  analysis  of  a  half*sine  acceleration  pulse 


ES 


(a)  Acceieration,  velocity,  and  displacement  time  histories 
Fig.  4  —  Output  for  analysis  of  FSP  acceleration  record 


(2^0  Oo«o)rQ)  Noiurol  Fr«qu«r^y  Hz 

(b)  Maximax  and  residual  shock  spectra 
Fig.  4  Output  for  analysis  of  FSP  acceleration  record 
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(d)  Tabular  listing  of  all  computed  data 
Fig.  4  Output  fur  analysis  of  FSP  acceleration  record 
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Appendix  1 

List  of  Program  Variables 


1  A  Arroy  for  rime  record  Oin(N0I 
2:A2  Zero  offset  •  Motrix  0in(7,2) 

30  Arroy  0rn(2) 

4  B3  Vorioble  used  in  subroutine  TIME 
5.04  Vorioole  used  m  suoroutme  TIME 
6 :C  Arroy  OIMtCBl 

7 : C0  Actual  number  of  analysts  frequencies  (defoul t  90) 

01  Cl  Element  of  E5  (f Irst  onolysis  f requency  per  frequency 
9tC2  Element  of  F5  (lost  onolysis  f requency  per  frequency 
10C3  Mox  possible  analysis  frequencies  (depends  on  memory 
110  I/O  device  oddress  fdefouli  32) 

12:01  -tie  number  on  exi,  tope  drive  to  store  do to 

1 3:03  File  number  on  ext.  tope  drive  to  read  do to  from 

14.E  Instontoneous  foe  tor  for  sample  step 
15;EI  Foctor  for  sample  step  slthin  frequency  ronge  ni-10Hz) 

16  E2  Foe tor  for  somple  step  within  frequency  ronge  n(10-100H2» 

17;E3  Foe  tor  for  somple  step  within  frequency  ronge  1 1 1 ( 1 00- 1 000H2 J 

^8.F  Instantaneous  volue  of  onolysis  frequency 

1 9 : F 1  Lowest  onol ys i s  f requency  (defoul t  1 ) 

20:F2  Highest  onolysis  frequency  (default  1000) 

21 :F3  ^owe^  onolysis  frequency  within  one  of  the  3  frequency  ronges 

22: f A  Upper  onolysis  frequency  within  one  of  the  3  frequency  ronges 

23:F5  Frequency  orroy  DIM(C0) 

24:F6  Lower  frequency  of  lin,  graph  (input) 

25!F7  Upper  frequency  of  lin.  groph  (input) 

26:F0  Cutoff  frequency  of  lowposs  filter 

27  G  Instontoneous  value  of  occeloration  f (2»p t »F) f2»X) 

28  G5  Accelerot’on  orroy  Din(C0) 

29  H  Time  increment  of  digitized  time  record 

30  HI  Exponent 

3'’H2  Time  increment  for  compuforion 
32 : I  Counter 

33  1 1  Counter  to  find  Cl 

34  12  Indicates  frequency  ronge  fl  to  III) 

35: 13  Somple  number  when  V2  occurs 

36.14  Somple  number  when  VI  occurs 
37  K  Target  or'^oy  for  groph  disp  lay 

38'AI  Foctor  of  F8  to  colculote  somple  frequency  within  frequency  ronqe  I 
e^oul t  2)  ^ 

39  A2  Focror  of  F8  to  calculate  somple  frequency  withm  frequency  range  (I 
def oul T  2) 

40  K3  Foctor  of  F8  to  colculote  somple  frequency  wlth;n  frequency  range  IT 

(def  oul t  1 0)  ^ 

41  L  Calibration  level  in  i npu  t  un i t  s 
42 . M  Number  of  doto  module 

43  ni  Mox t mum  of  orroy  \  for  outoscole 

44  h9  ColibroTion  Moirix  0IMf7,21 

45  N0  Number  of  doto  points  on  record  (defoult  4096) 

46:N  Number  of  doto  points  considered  dur i ng  colculotron 

47  Nl  _ower  veritcol  level  of  log  groph  disploy 

48  N2  Upcer  very. col  level  of  log  groph  disploy 

49. N4  Number  of  doto  points  at  the  beginning  of  me  residuol  i-esporse 
50tN5  Number  of  doto  points  after  rrunctot.on 
5’  N6  Lower  frequency  of  log  groph  disploy 
52  N7  Upper  frequency  of  log  groph  display 

53. C  Indicotes  me  port  of  the  progrom  currently  looded  (i  or  2) 

54P  Selects  g- . d  on  log.  disploy  (P>I0  no  grid,  defoult  10) 


ronge) 
ronge) 
ov  a i 1 ob  t  e ) 


25 


55. PI  Error  iodicaror  if  >0  902 

5B  Q  Number  of  periods  considered  per  analysis  frequency 

57  K  Selecis  riAXIMAX-  ond/or  RESIDUAL  SHOCK  SPECTRun ( Keys  6  to  9l 

58  R I  Selecis  d»v.  linear  spectre  (user  defined  Keys  16  to  23) 

59  R2  Ind  cores  d » splocemen t ( J ) , veloc i t y (2)  or  accelerot ion(3)  record 

60  R3  Equcls  X3  before  N4  (initioll 
6-  R4  Equcls  X3  of  ter  Iresidual) 

62  R5  Flog  for  offset -procedure 

63 :R6  Flog  for  deleting  colifcrotion  routjre 
54  R8  Flog  >n  col'brotion  routine 

65  S  Colibrotion  factor 

66  SI  Colibrotion  factor  after  a i f f eren r i o t i on  ond/or  integrotion 
67:S2  Coi-broiion  foctor  of  ter  E-procedure 

63  T  I  ns T on t oneous  volue  of  time  of  maximum  (T3  or  T4) 

69-71  Beginning  of  time  record 

70  :T3  Time  of  mox .  response  before  N4  UnitioJ) 

7'  t4  Time  o^  mox.  response  Of  ter  N4  (residual) 

72  T5  Time  orroy  (conto*ns  max  volues  of  either  T3  or  T41  Orri(C0) 

73  U  Time  difference  m  decimal  value 

74 . V  Instontoneous  volue  of  pseudvel oc i t y  ({2*pifF)*X) 

75 : V •  hox i mum  of  orroy  A 
76  V2  f1  i  n  I  mum  of  orroy  A 
77:V5  Velocity  orroy  Din(C0) 

78  ^  Foctor 

79  w  Number  of  spectra  pro  plot  (flul  t  i  Display  node)  and  flag  for  ws-qisplo 

y 

80  W2  Array.  Din  W2(WI ) 

8):X  X-coordino:e  of  pointer 

82  X  Instontoneous  value  of  displocement  (maximum  of  R3  or  R4 » 

83:X’  Deflection  ot  -0.5»h2 
84:X2  Deflection  ot 

65:  X3  riox .  response  of  onolysis  frequency  P 

86:X5  Displacement  orroy 

87 ; Y  ^-coordinote  of  pointer 

88  Y  ! ns r on r oneous  value  of  of  mpx .  response  velocity  of  ter  N4 

8S  r  y5  Velocity  orroy  DIH (C01 

90 :Z  I rs ton T oneous  value  of  of  phase 

9l;Z5  Phase  orroy  DIMfCB) 

92:AS  String  for  groph  dispioy 

93:3*  System  time  when  process -procedure  i s  comp  1 e t ed 
94 : C*  Time  S  t  r I ng 

95 . D*  String  :o  address  data  module 
96.3*  String  for  g^oph  dispioy 
97 . H*  String  for  groph  d  <  sp 1 oy 
98  L*  String  ^0  check  cond.’ion  of  data  module 
99 • n*  String  for  groph  dispioy 
1 00  N*  String  for  groph  dispioy 
I  0 1  . P*  Spring  for  gr oph  display 
102  Q*  Spring  C  f  *0"  stops  onolysis  procedure) 

03 :R*  “NRL  -  CODE  5842“ 

104:T*  Contains  system  time 
05:V*  Contoins  text  for  vert  level 
'  06  W*  Title  to  identify  record 
I07;X*  Segment  string 
'09  Y$  Segment  strtog 
!09  Z*  Segment  string 


Appendix  2 
Program  Listing 

1  REM  ////////////////////////  eSPRGM/PART  I  \\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2  INIT 

3  GO  TO  260 

4  IF  0=1  THEN  92 

5  GO  TO  7450 

8  IF  0=1  THEN  92 

9  GO  TO  7530 

12  IF  0=1  THEN  92 

13  CO  TO  8980 
16  GO  TO  190 
20  GO  TO  100 

24  IF  0=1  THEN  92 

25  R»1 

26  GO  TO  7690 

28  IF  0=1  THEN  92 

29  R»2 

30  GO  TO  7690 

32  IF  0=1  THEN  92 

33  R=3 

34  GO  TO  7690 

36  IF  0=1  THEN  92 

37  GO  TO  9510 

40  IF  0=1  THEN  92 

41  GO  TO  9780 

44  IF  0=2  THEN  90 

45  CHARSIZE  3 

46  IF  R6=0  THEN  94 

47  IF  R6=l  THEN  4200 

48  IF  0=2  THEN  90 

49  GO  TO  530 

50  DELETE  A 

52  IF  0=2  THEN  90 

53  GO  TO  2510 

54  DELETE  A 

56  IF  0=2  THEN  90 

57  GO  TO  5250 

60  IF  0=2  THEN  90 

61  R5=l 

62  CHARSIZE  3 

63  GO  TO  2410 

64  IF  0=1  THEN  92 

65  Rl=l 

66  GO  TO  6010 

68  IF  0=1  THEN  92 

69  Rl=4 

70  GO  TO  6010 

72  IF  0=1  THEN  92 

73  Rl=2 

74  GO  TO  6010 

76  IF  0=1  THEN  92 

77  Rl*3 

78  GO  TO  6010 

80  IF  0=1  THEN  92 

81  Rl=5 

82  GO  TO  6010 

90  PRINT  "You  got  to  PRESS  KEY  5  first !G” 

91  RETURN 

92  PRINT  "You  got  to  PRESS  KEY  4  firstIG" 

93  RETURN 

94  PRI  "LJJTHE  CALIBRATION  PROCEDURE  HAS  BEEN  DELETED  IN  ORDER  TO  SAUE"; 
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95  PRINT  “  MEMORY.  IF  YOU  WANT  TO  REDEFINE  THE  CALIBRATION  MATRIX  M9  ”  ; 

96  PRINT  "RELOAD  THE  PROGRAM  (PRESS  USER  DEFINED  KEY  4  AND  5  IN  A  ROW* 

97  PRINT  "OR  PRESS  >  AUTOLOADO  *  GG  " 

99  RETURN 

10  0  REM  LOAD  ROUTINE  xxxxxxxxx*x**xx*x**»***x»x»**«x*»»x 

no  IF  0»l  THEN  160 

120  DELETE  F5 , X5 > U5 , G5 , T5 , Y5 . Z5 , C , K  .  A 
130  R6=l 
140  0=1 

150  APPEND  "eSPRGM/PARTI’*;l 

160  PRINT  "LGSPRGM/PART  I  IS  LOADED!" 

170  DELETE  5380,20000 
180  RETURN 
190  IF  0=2  THEN  240 
200  0=2 

210  DELETE  270.20000 

220  APPEND  '*eSPRGM/PARTII";260 

230  RENUMBER  6000,10,260 

240  PRINT  "LGSPRGM/PART  II  IS  LOADED!" 

250  RETURN 

260  REM  X  X  xxxxxxxxx  xxxx  DEFAULT  CONDITION  xxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 
270  W$="TEST" 

280  0=1 
290  W1=0 
300  R5=0 
310  R6=l 
320  D=32 
330  Fl=l 
340  F2=1000 
350  F8=1000 
360  Nl=0.01 
370  N2=100 
380  N6=l 
390  N7=1000 
400  C0=90 
410  N0=4096 
420  Kl=2 
430  K2=2 
440  K3=10 
450  P=10 

460  DIM  M9(7,2) ,C(C0 ) ,A2<7,2) ,R$( 15) ,T$( 18) ,B$( 18) ,C$< 18) 

47  0  4>IM  Di(8)  ,H$(20)  .yi(22)  ,G$(7)  ,M$(7)  ,N$<8)  ,X$(22)  ,PiK35)  ,L$<4)  ,Q$<  1  ) 
480  R^>="NRL  -  CODE  5842" 

490  M9^ 0 
500  A2=0 

510  PRINT  "  GSPRGM/PART  I  IS  LOADED  (DEFAULT  CONDITIONS)!" 

520  RETURN 

530  REM  X xxxx« xxxx-xxxxxxxx  READ  RECORD  xx-xxxxxx-xxxxxxxxxxxxxxxxxxxxx-xxxx 
540  DELETE  4200,5240 
550  R6=0 

560  PRINT  eD,17: 1 .2, 1 .3 
570  CHARSI2E  3 

580  PRI  "LTHE  CURRENT  TITLE  IS  J  ENTER  A  NEW  TITLE  IF  NECESSA"; 

590  PRINT  "RY.  OTHERWISE  PRESS  'RETURN'. JG" 

600  INPUT  A$ 

610  IF  Ai=""  THEN  630 
620  W4=AS 
630  DIM  A(N0) 

640  CHARSI2E  3 

650  IF  R5=0  THEN  660 
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660  U3«130 

670  IF  D<  >4  THEN  690 
680  U3»150 

690  PRINT  "JENTER  NUMBER  OF  DATA  MODULE  TO  READ  FROM:G"; 

700  INPUT  M 

710  IF  R5*l  THEN  740 

720  IF  M9<M,2)==0  THEN  2390 

730  S=M9(M,2) 

740  PRINT  86:‘’P" 

750  GO  TO  M  OF  760,780,800,820,840,860,880 
760  D$«"I01B4096'' 

770  CO  TO  890 
780  D*»“I02B4096‘* 

790  GO  TO  890 
800  D$=”I03B4096” 

810  GO  TO  890 
820  D^*’*I04B4096‘' 

830  GO  TO  890 
840  DS="I05B4096” 

850  GO  TO  890 
860  D»="I06B4096” 

870  GO  TO  890 
880  D$="I07B4096” 

890  PRINT  e6:Di 
900  INPUT  e6:A 
910  P^s"D'* 

920  PRINT  e6:P$ 

930  INPUT  ^6i?% 

940  AaA-A2(M,2) 

950  HS="SA’* 

960  PRINT  96iH% 

970  INPUT  e6;H$ 

980  H1=UAL(H*) 

990  IF  HI >12  THEN  1180* 

1000  GO  TO  HI  OF  1,1,1,1010,1030,1050,1070,1090,1110,1130,1150,1170 

1010  HsO.005 

1020  GO  TO  1200 

1030  H=0.01 

1040  GO  TO  1200 

1050  H*0.025 

1060  GO  TO  1200 

1070  H=0.05 

1080  GO  TO  1200 

1090  H=0 . 1 

1100  GO  TO  1200 

1110  H=0.25 

1120  CO  TO  1200 

1130  H«0.5 

1140  GO  TO  1200 

1150  H*1 

1160  GO  TO  1200 

1170  H*2.5 

1180  PRINT  "ENTER  TIME  INCREMENT  IN  MS 

1190  INPUT  H 

1200  IF  R5al  THEN  2430 

1210  CHARSIZE  2 

1220  CALL  "TIME^T^ 

1230  PRINT  6D: "LK" jR4 
1240  PRINT  8D;Ti 
1250  UINDOU  0  ,U3,0 ,100 
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1260  yiEUPORT  0,U3, 0,100 
1270  IF  W1 >0  THEN  1350 
12B0  PRINT  eD,17!2,3 
1290  CHARSIZE  4 
1300  MOVE  8D:W3/2,98 
1310  FOR  1=1  TO  LEN(U*)/2 
1320  PRINT  @0: "H" ; 

1330  NEXT  I 
1340  PRINT  eD;W* 

1350  PRINT  eD,17!l .2,1 .8 

1360  CHARSIZE  2 

1370  Tl=H»AEiS(VAL<P*) ) 

1300  GO  TO  R2  OF  1390,1520,1660 

1390  GOSUEi  1840 

1400  S1=S 

1410  GOSUB  2020 

1420  PRINT 

1430  CALL  "DIF3",A,A 

1440  S1=S1/H 

1450  GOSUB  1900 

1460  GOSUB  2020 

1470  CALL  "DIF3",A,A 

1480  S1=S1 »10*3/(H*9.8044) 

1490  GOSUB  1960 
1500  GOSUB  2020 
1510  GO  TO  1780 
1520  GOSUB  1900 
1530  S1=S 
1540  GOSUB  2020 
1550  PRINT 

1560  CALL  "DIF3",A,A 
1570  S1=S1»10*3/(9.8044*H) 

1580  GOSUB  1960 

1590  GOSUB  2020 

1600  CALL  "INT  ”,A,A 

1610  CALL  "INT  ",A,A 

1620  S1=S1»1 . 0E-3»H*2»9.8044 

1630  GOSUB  1840 

1640  GOSUB  2020 

1650  GO  TO  1780 

1660  GOSUB  1960 

1670  S1=S 

1680  GOSUB  2020 

1690  PRINT  8D."*" 

1700  CALL  "INT", A, A 

1710  S1=S1»1 .0E-3»H»9.8044 

1720  GOSUB  1900 

1730  GOSUB  2020 

1740  CALL  "INT", A, A 

1750  S1=S1*H 

1760  COSUB  1840 

1770  GOSUB  2020 

1780  VIEWPORT  0,U3, 0,100 

1790  WINDOW  0,W3, 0,100 

1800  MOVE  eD;0,0 

1810  PRINT  @D; "*  oriq.  tiwe  record 

1820  delete  a 

1830  RETURN 

1840  VIEWPORT  5,W3-30,4,32 
1850  V** "DISPLACEMENT" 
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1860  nax  : 

1870  nin 

1880  N$=5“  C«rt]" 

1390  RETURN 

1900  VIEWPORT  5,W3-30,36,64 
1910  V^=”VELOCITY" 

1920  G*=s"V  «ax  :  '• 

1930  H$="V  Min  ; 

1940  N$='*  [M/sec]  ” 

1950  RETURN 

1960  VIEWPORT  5  ,  W3-'30 , 68 , 96 
1970  V$=''ACCELERATION" 

1980  GS=‘'A  Max  ;  ” 

1990  M5=‘’A  Min  ;  ” 

2000  N$=''  Cgl" 

2010  RETURN 

2020  CALL  "MAX”, A, VI ,14 
2030  CALL  "MIN", A, V2, 13 
2040  WINDOW  0,N0,V2,V1 
2050  IF  D=32  THEN  2110 
2060  MOVE  60:0,0 
2070  FOR  1=1  TO  NO 
2080  DRAW  6D:I,A<I) 

2090  NEXT  I 

2100  GO  TO  2120 

2110  CALL  "DISP",A 

2120  AXIS  60:400,0 

2130  IF  R5=l  THEN  2340 

2140  FOR  1=0  TO  NO  STEP  2000 

2150  MOVE  6D;I,0 

2160  PRINT  60; "HJ",I»H; 

2170  NEXT  I 

2180  PRINT  60: "HHHHHHHHHHHCMS3 " 

2190  MOVE  eO:4200,vT 

2200  PRINT  60;  USING  2240 : "J" ,G$,V1 xSl ,N$; 

2210  MOVE  60;4200,V1 
2220  PRINT  60; "J"; 

2230  PRINT  60;  USING  2240  :  "  J  "  ,  M*  ,  V2^tSl  ,  N$  > 

2240  IMAGE  A » 7A » 60 . 1 0 » 1 1 A 
2250  MOVE  60: 0 ,0 .5x(V2+Vl ) 

2260  PRINT  60: "HH"  ; 

2270  FOR  1=1  TO  LEN<V$)/2 
2280  PRINT  60; "K"; 

2290  NEXT  I 

2300  FOR  1=1  TO  LEN<V$) 

2310  X$=SEG(V$, 1 , 1  ) 

2320  PRINT  60;X*;"HJ"; 

2330  NEXT  I 
2340  RETURN 

2350  PRINT  6D;"JDATA  MODULE  ")Mj"  IS  IN  GROUND  MODE  MGG" 
2360  RETURN 

2370  PRINT  "JDATA  MODULE  ";M;"  IS  IN  CALIBRATION  MODE  !!GG" 
2380  RETURN 

2390  PRINT  "JREAD  THE  CALIBRATION  STEP  FIRST !C" 

2400  RETURN 

2410  PRINT  "LOFFSET  -  CORRECTION;" 

2420  GO  TO  630 
2430  VIEWPORT  0,110,10,90 
2440  CALL  "INT",A,A 
2450  GOSUB  2020 
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2460  POINTER  X>Y,a$ 

2470  A2(M,1 )«M 

2480  A2<M,2)«Y/X+A2<M>2) 

2490  R5«0 
2500  RETURN 

2510  REM  PROCESS 

2520  CMARSIZE  3 
2530  CALL  “TIME’^BS 
2540  PRINT  @0:  “L'*jB$)'*J" 

2550  PRINT  e32,26;2 
2560  IF  F1<1  THEN  2620 
2570  IF  FI >1000  THEN  2620 
2580  IF  F2<1  THEN  2630 
2590  IF  F2>1000  THEN  2630 
2600  IF  F2<=F1  THEN  2630 
2610  GO  TO  2660 
2620  IF  C0=*1  THEN  2660 

2630  PRI  “ANALYSIS  FREQUENCY  IS  OUT  OF  RANGE  ENTER  NEW  VALUES  OF  F1,F2G" 
2640  INPUT  F1,F2 
2650  GO  TO  2530 

2660  DELETE  F5 , X5 , V5 ,G5 , T5 , Y5 > Z5 , C > K , A 
2670  DIM  A(N0) 

2680  C3aINT<<MEM0RY/8-l8)/8) 

2690  DIM  F5<C0),X5(C0>,V5<C0>,G5(C0),T5(C0),Y5(C0>,25(C0),C<C0) 

2700  ElaINT<  1000/<KlitF8*H)) 

2710  E2*INT<  1000/(K2iiF8itH)  ) 

2720  E3*INT(  1000/(K3^^F8«H)  > 

2730  F5<1>=F1 

2740  IF  C0=1  THEN  2780 

2750  FOR  1=2  TO  CO 

2760  F5(I)=<F2/F1 >^(1/(C0-1) )«F5(I-1 > 

2770  NEXT  I 
2780  X5»0.01 
2790  V5=0.01 
2800  G5=0.01 
2810  T5=0.01 
2820  Y5=0.01 
2830  Z5=0.01 
2840  Il«l 
2850  12=0 

2860  IF  FI <10  THEN  2900 
2870  12=1 

2880  IF  FUlOO  THEN  2900 

2890  12=2 

2900  GO  TO  3070 

2910  E*E1 

2920  GOSUB  3160 

2930  IF  II >C0  THEN  3010 

2940  GO  TO  3070 

2950  E»E2 

2960  GOSUB  3160 

2970  IF  I1>C0  THEN  3010 

2980  GO  TO  3070 

2990  E*E3 

3000  GOSUB  3160 

3010  DELETE  A 

3020  CALL  “TIME‘\T$ 

3030  GOSUB  4030 
3040  PRINT  032,26:0 

3050  PRINT  0D:“  MAX.  POSSIBLE  ANALYSIS  FREQUENCIES;  *';C3 


3060  RETURN 
3070  F3=F5(I1) 

3080  Cl=Il 

3090  I2»I2+1 

3100  FOR  11*1  TO  CO 

3110  IF  F5<I1)>10^I2  THEN  3130 

3120  NEXT  II 

3130  F4*F5(I1-1) 

3140  GO  TO  12  OF  2910,2950,2990 

3150  RETURN 

3160  PRINT  86: ”P”;D$ 

3170  INPUT  96; A 
3180  A*A-A2(M,2) 

3190  IF  E<=1  THEN  3260 
3200  J=1 

3210  FOR  1=1  TO  NO  STEP  E 
3220  A(J>=A(I> 

3230  J*J+1 
3240  NEXT  I 
3250  GO  TO  3270 
3260  E=1 
3270  N5=N0/E 
3280  H2*H*E 

3290  GO  TO  R2  OF  3300,3340,3370 

3300  CALL  ”DIF3”,A,A 

3310  CALL  "DIF3",A,A 

3320  S2*S 

3330  GO  TO  3380 

3340  CALL  ”DIF3”,A,A 

3350  S2*S*H2 

3360  GO  TO  3380 

3370  S2*S*10^-3»H2^2*9.8044 

3300  C2=I1-1 

3390  P1»0.0062831853^^F3^^H^ 

3400  FOR  I*C1  TO  C2 
3410  C(I)*2»C0S(P1 ) 

3420  IF  C1*C2  THEN  3440 
3430  P1=P1»(F4/F3>^(1/(C2-C1 )  ) 

3440  NEXT  I 

3450  IF  Pl<0.902  THEN  3500 
3460  IF  143.566/H2>F2  THEN  3500 

3470  PRINT  "LGUARNING;  FREQUENCIES  ARE  TOO  HIGH  FOR  YOUR  TIHE  INCREMENT" 
3480  PRINT  "  ERRORS  ARE  OUTSIDE  THE  RANGE  -7  TO  ^3  PERCENT  FOR" 

3490  PRINT  "  FREQUENCIES  ABOVE  *M  43 . 566/H2 ; "  HERTZ" 

3500  FOR  I=C1  TO  C2 
3510  F»F5<I) 

3520  IF  F<=1  THEN  3550 
3530  N=1  000/(H2»5<*LGT(F)  ) 

3540  IF  N<N5  THEN  3560 

3550  N=N5 

3560  X1*0 

3570  X2»0 

3580  R3*0 

3590  R4-0 

3600  IF  F*0  THEN  3630 
3610  N4*N-l-“500/<F»H2> 

3620  IF  N4>3  THEN  3640 
3630  N4*3 

3640  FOR  J*1  TO  N-1 
3650  X3»C( I)*X2-X1>A( J) 
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3660  IF  J<N4  THEN  3710 

3670  IF  ABS(X3)<ABS<R4)  THEN  3740 

3680  R4=X3 

3690  T4»J 

3700  GO  TO  3740 

3710  IF  ABS<X3)<ABS(R3)  THEN  3740 

3720  R3=:X3 

3730  T3=J 

3740  X1=X2 

3750  X2=X3 

3760  NEXT  J 

3770  IF  ABS(R3) >A&S(R4)  THEN  3800 

3780  R3=R4 

3790  T3=T4 

3800  X=S2*R3 

3810  y=0 , 00628318»F)frX 

3820  G=*6.2883ia*F*V/9.8044 

3830  T=Tl+(T3-0 .5)»H2 

3840  Y=0.006283l8itF)tS2xR4 

3850  Z=1  .  0E-3it(Tl  +  (T4-0 .5)  )eH2)»F  +  0 .25 

3860  Z=:Z-INT(Z) 

3870  IF  Y>0  THEN  3890 
3880  Z=Z-0.5 

3890  IF  Z<=0.5  THEN  3910 

3900  Z=Z-'l 

3910  Z=360)*Z 

3920  F5<I)=F 

3930  X5(I)=X 

3940  y5(I)=y 

3950  G5(I)=G 

3960  T5(I)=T 

3970  Y5(I)=Y 

3980  Z5(I)=Z 

3990  PRINT  80:  USING  4000 : I ,F,U,N; 

4000  IhAGE  20T,3D,4X,4D. 1D,4X,5D.2D,3X,6D 
4010  NEXT  I 
4020  RETURN 

4030  REM  SUBROUTINE  TIME 

4040  DIM  B(2) 

4050 

4060  FOR  1*1  TO  2 
4070  X^*SEG(Ci, 1 1 ,8) 

4080  Y$*SEG(XS,4,2) 

4090  Z$*SEG(X$,7,2) 

410  0  B(  I  )=UAL(Xi)  +  (VAL(  Y»)  »60-*-VAL(Z$)  )/360  0 
4110  C$*T$ 

4120  NEXT  I 
4130  U*B(2)“B( 1 ) 

4140  B3*<U~INT<U) )*60 
4150  B4*(B3-INT(B3))»60 

4160  PRINT  eD:  ••^CALCULATION  TIME  (hr  ;  Min  :  sec  )  ;  "  ; 

4170  PRINT  80;  USING  4180  :  INT(U>  •MNT<B3)  , '•:*•,  INT<B4) 

4180  IMAGE  2D, A, 2D, A, 2D 
4190  RETURN 

4200  REM  xxxxxxxxxxxxxxx  CALIBRATION  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
4210  DIM  A(N0) 

4220  CHARSIZE  3 

4230  PRINT  "LARE  YOU  DEALING  WITH  J  A  DISPLACEMENT  RECORD  (MM>"; 

4240  PRINT  -  (ENTER  1)  A  VELOCITY  RECORD  (M/SEC>  (ENTER  2)  OR  A‘*j 

4250  PRINT  "  ACCELERATION  RECORD  (G)  (ENTER  3)5  •• ; 


4260  INPUT  R2 

4270  PRI  •• _ DO  YOU  WANT_^  TO  READ  A  CALIBRATION  STEP  (ENTER  1)  OR  “ ; 

4280  PRINT  "  TO  ENTER  CALIBRATION  FACTOR  (ENTER  2>G 
4290  INPUT  R8 
4300  M9=0 

4310  IF  R8=2  THEN  4960 
4320  FOR  rt=l  TO  7 

4330  PRI  "LENTER  CALIBRATION  LEVEL  FOR  CHANNEL  ”  ; M ; ”  IN  INPUT  UNIT5G:"; 

4340  INPUT  L 

4350  GOSUB  4750 

4360  IF  S4=0  THEN  4390 

4370  IF  S4=7  THEN  4420 

4380  GO  TO  4450 

4390  GOSUB  2310 

4400  CALL  "WAIT", 2 

4410  GO  TO  4730 

4420  GOSUB  2330 

4430  CALL  "WAIT”, 2 

4440  GO  TO  4730 

4450  PRINT  @6; "P" 

4460  GO  TO  M  OF  4470 , 4490 > 45 1  0 , 4530 > 4550 , 4570 , 4590 
4470  Df='*I01B4096" 

4480  GO  TO  4600 
4490  D$=”I02B4096" 

4500  GO  TO  4600 
4510  D$="I03B4096" 

4520  GO  TO  4600 
4530  D$='*I04B4096” 

4540  GO  TO  4600 
4550  D$= ”10584096” 

4560  GO  TO  4600 
4570  Di=”I06B4096" 

4580  GO  TO  4600 
4590  D^*”I07B4096” 

4600  PRINT  e6:D$ 

4610  INPUT  @6: A 

4620  WINDOW  0,5000,-512,512 

4630  VIEWPORT  0,130,0,100 

4640  PRINT  "LCHANNEL  ”;M;”  JJSET  CROSSHAIR  TO  UPPER  CALIBRATION  LEVELG" 
4650  CALL  ”DfSP”,A 
4660  AXIS 

4670  POINTER  XI  ,Y1  ,Q$ 

4680  PRINT  ”*TJJ  AND  LOWER  CALIBRATION  LEVELG !  ” 

4690  POINTER  X27Y2,Q$ 

4700  IF  Y2>Y1  THEN  5230 

471  0  H9<M,  1  )=»rt 

4720  M9(M,2)*L/(Y1-Y2> 

4730  NEXT  M 
4740  RETURN 

4750  GO  TO  M  OF  4760,4780,4800,4820,4840,4860,4880 
4760  L$=”A01" 

4770  GO  TO  4890 
4780  L$*”A02” 

4790  GO  TO  4890 
4800  Li=”A03” 

4810  GO  TO  4890 
4820  L$-”A04” 

4830  GO  TO  4890 
4840  Li=”A05” 

4850  GO  TO  4890 
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4860  L$=“A06'* 

4870  GO  TO  4890 
4880  Li="A07“ 

4890  PRINT  e6;Li 
4900  INPUT  @6;L^ 

4910  Sl=INT<VAL<Li)/l 000) 

4920  S2=INT< (UAL(L$)-S1»1000)/100) 

4930  S3=INT( (UAL(L$)-S1*1000"S2*100 VIO) 

4940  S4=UAL ( L$) -SI  *1  00 0--S2»l  0  0-53*10 

4950  RETURN 

4960  PRINT  "LENTER;" 

4970  FOR  h=l  TO  7 

4980  GOSUB  4750 

4990  IF  S4=0  THEN  5020 

5000  IF  S4=7  THEN  5040 

5010  GO  TO  5060 

5020  GOSUB  2310 

5030  GO  TO  5200 

5040  GOSUB  2330 

5050  GO  TO  5200 

5060  PRI  "JCALIBRATION  FACTOR  ( INPUT-UNITS/ 1 UOLT )  FOR  CHANNEL 
5070  INPUT  L 

5080  GO  TO  S4  OF  5090 , 5 1 1 0 , 51 30 > 5 I 50 , 5 1 70 , 51 90 

5090  M9(M,2)=L/512 

5100  GO  TO  5200 

5110  M9(M,2)=L*0 .2/512 

5120  GO  TO  5200 

5130  M9(M,2)=L*5/512 

5140  GO  TO  5200 

5150  M9(M,2)=L*0 . 1/512 

5160  GO  TO  5200 

5170  M9(M,2)=L*2/512 

5180  GO  TO  5200 

5190  M9(M,2)=L*0 .5/512 

5200  M9(M>1>=M 

5210  NEXT  M 

5220  RETURN 

5230  PRINT  ""JJJJJ  UPPER  CALIBRATION  LEVEL  FIRST!!  TRY  IT  AGAIN! !G" 

5240  GO  TO  4670 

5250  REM  »»»**»* *****x*x»»x****)t  UR  ITE  RECORD  ************************** 
5260  DIM  A(N0) 

5270  CHARSIZE  3 

5280  PRINT  "LENTER  FILE  NUMBER  (DEVICE  NR. 3)  TO  WRITE  RECORD  ONG s 
5290  INPUT  D3 
5300  FIND  83:D3 
5310  PRINT  @6:  ‘»P'*;D$ 

5320  INPUT  e6:A 
5330  A=A-A2(M,2) 

5340  WRITE  83 : W^ > R2 , H , T 1 , M , M9 ( M , 2 ) , A2 ( M » 2 ) , A 
5350  PRINT  "JRECORD  IS  STORED  IN  FILE  *0  ;  D3 
5360  DELETE  A 
5370  RETURN 


6000  REM  SPRGM/PART  II  FILE  3 

6010  REM  *************** ********  LIN .  SPECTRA  ************************** 
6020  DELETE  K 
6030  DIM  K<C0) 

6040  At=:"MAXIMAX  SHOCK  SPECTRUM" 

6050  GO  TO  R1  OF  6060 > 6080 > 61 00 , 61 20 , 61 50 
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6060  K=ABS(V5> 

6070  GO  TO  6170 
6080  K»ABS(G5) 

6090  GO  TO  6170 
6100  K»ABS(X5> 

6110  GO  TO  6170 

6120  Af»^'RESIDUAL  SHOCK  SPECTRUM” 

6130  K»ABS<Y5> 

6140  GO  TO  6170 
6150  A4»”PHASE” 

6160  K^ZS 

6170  CALL  "MAX”, K, Ml ,14 
6180  CHARSIZE  3 

6190  PRINT  "LJFREQUENCY  RANGEG: 

6200  INPUT  F6,F7 
6210  IF  D=4  THEN  6240 

6220  VIEWPORT  8*1 . 79 , 1 30-4^t  1 . 79 , 34^2 . 82 , 1  0  0-3it2 . 82 
6230  GO  TO  6260 

6240  VIEWPORT  8*1.79,150-4*1.79,3*2.82,100-3*2.82 
6250  PRINT  eD,17; 1 .2,1 .8 
6260  CALL  ••TIME",T$ 

6270  CHARSIZE  1 
6280  PRINT  §D; ”L";R$ 

6290  PRINT  eD:T$ 

6300  IF  Rl=5  THEN  6510 

6310  IF  ABS(M1)<10  THEN  6350 

6320  IF  ABS<M1)<100  THEN  6390 

6330  IF  ABS(hl)<1000  THEN  6430 

6340  IF  ABS(M1>>1000  THEN  6470 

6350  M1*INT(M1 )+l 

6360  WINDOW  F6,F7,0,M1 

6370  AXIS  8D: <F7-F6>/10,0 .5 

6380  GO  TO  6580 

6390  M1=INT(M1/10)*10+10 

6400  WINDOW  F6,F7,0,H1 

6410  AXIS  eD: (F7-F6)/10,5 

6420  GO  TO  6580 

6430  M1=INT(M1/100)*100+100 

6440  WINDOW  F6,F7,0,M1 

6450  AXIS  8D: (F7-F6)/10,50 

6460  GO  TO  6580 

6470  M1*INT<M1/1000)*1000+1000 

6480  WINDOW  F6,F7,0,M1 

6490  AXIS  eD: (F7-F6>/10,500 

6500  GO  TO  6580 

6510  WINDOW  F6,F7, 0,180 

6520  AXIS  eD: (F7-F6)/10,45 

6530  HOVE  8D:F6,180 

6540  DRAW  8D:F7,180 

6550  DRAW  eD:F7,0 

6560  MOVE  CD; (F7+F6)/2,180 

6570  CO  TO  6610 

6560  HOVE  9D:F6,M1 

6590  DRAW  tD:F7,Ml 

6600  DRAW  9D!F7,0 

6610  PRINT  9D, 17;2,3 

6620  CHARSIZE  4 

6630  GOSUB  6700 

6640  IF  W1>0  THEN  6800 

6650  A9«W$ 
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6660  PRINT  eD>17; I .2,1 .8 

6670  CHARSIZE  2 

6680  GOSUB  6700 

6690  GO  TO  6800 

6700  IF  Rl=5  THEN  6730 

6710  MOVE  @D:  <F7'^F6)/2,M1 

6720  GO  rO  6740 

6730  eO;  (F7>F6)/2, 180 

6740  PRINT  @D: "KK"  > 

6750  FOR  1=1  TO  LEN(A$)/2 
6760  PRINT  @D:  ; 

6770  NEXT  I 
6780  PRINT  eD;A$; 

6790  RETURN 

6800  PRINT  eD,17: 1 .2,2 
6810  CHARSIZE  2 

6820  H$='*Natural  Frequency  Hz” 

6830  MOVE  SD; (F7+F6)/2,0 
6840  PRINT  eD: "JJJ”  j 
6850  FOR  1=1  TO  LEN(H$)/2 
6860  PRINT  eD; ”H"  j 
6870  NEXT  I 
6880  PRINT  eD:H^; 

6890  GO  TO  R1  OF  6900,6940,6960,6920,6980 
6900  U$=”Pseudovelocity  M/sec” 

6910  GO  TO  7010 

6920  U$=”Rel.  Velocity  M/sec" 

6930  GO  TO  7010 

6940  U$=”Ab5.  Acceleration  g '' 

6950  GO  TO  7010 

6960  U$=”Rel.  Displacewent  mm" 

6970  GO  TO  7010 
6980  Phase  deg” 

6990  HOVE  eD:F6,90 

7000  GO  TO  7020 

7010  MOVE  eD:F6,0 .5*INT(M1 ) 

7020  PRINT  eD;"HHHHHH"; 

7030  FOR  1=1  TO  LEN<V$)/2 
7040  PRINT  eD; ”K”; 

7050  NEXT  I 

7060  FOR  1=1  TO  LEN(V$) 

7070  X$=SEG(V$,I,1  ) 

7080  PRINT  eD;X$; ”HJ”  j 
7090  NEXT  I 

7100  FOR  I=F6  TO  F7  STEP  <F7-F6)/10 

7110  MOVE  eD:I,0 

7120  PRINT  eD; "HJ”;INT(I); 

7130  NEXT  I 

7140  IF  Rl=5  THEN  7270 
7150  IF  ABS(M1)<=10  THEN  7190 
7160  IF  ABS(M1><=100  THEN  7210 
7170  IF  ABS(M1)<=1000  THEN  7230 
7180  IF  ABS<M1)>1000  THEN  7250 
7190  FOR  1=0  TO  ABS(Ml)  STEP  1 
7200  GO  TO  7280 

7210  FOR  1=0  TO  ABS(Ml)  STEP  10 
7220  GO  TO  7280 

7230  FOR  1=0  TO  ABS(Ml)  STEP  100 
7240  GO  TO  7280 

7250  FOR  1=0  TO  ABS<M1)  STEP  1000 


GO  TO  7280 

FOR  1*0  TO  180  STEP  45 
MOVE  eD5F6,I 
PRINT  SD;  “HHH'*  ; 

NEXT  I 

MOVE  eD:F6;0 
PRINT  90,17:1 ,2 

(Zero  Damping >  ’* 

PRINT  90: 

1*1 

FOR  I*I  TO  CO 
IF  F5(I)>0.01  THEN  7390 
NEXT  I 

MOVE  eOiFSd)  ,ABS<K(I)  ) 

FOR  I=I  TO  CO 
IF  F5(I)*0.01  THEN  7360 
DRAW  90:F5<I) ,ABS(K(I) ) 

NEXT  I 
RETURN 

REM  WRITE  OAT  A  «««««««««««««««« > 

CHARSIZE  3 

PRINT  "JLGENTER  FILE  NUMBER  (DEVICE  NR. 3)  TO  WRITE  DATA  ON: 

INPUT  01 
FIND  93:01 

WRITE  93 : CO , W$ , F5 , X5 , V5 , G5 , T5 , Y5 , Z5 
PRINT  "RESULTS  STORED  IN  FILEG  ♦";D1 
RETURN 

REM  READ  DATA 

CHARSIZE  3 

PRINT  "JLGENTER  FILE  NUMBER  (DEVICE  NR. 3)  TO  READ  DATA  FROM:  "> 
INPUT  03 

DELETE  F5,X5,V5,G5,T5,Y5,25,K,A 
FIND  93;D3 
READ  93: CO 
FIND  93:03 

DIM  F5(C0) ,X5(C0) ,V5(C0) ,G5(C0) ,T5(C0) ,Y5(C0) ,Z5(C0) ,K(CO) 

READ  93 : CO , W$ , F5 , X5 , V5 , G5 , T5 , Y5 , Z5 

IF  W1 >0  THEN  7650 

PRINT  "DATA  ARE  IN  MEMORY!G“ 

RETURN 

REM  LOG.  SPECTRA 

VIEWPORT  5*1 .79,110-4*1 .79,10+3*2.82,100-3*2.82 
GO  TO  7740 

VIEWPORT  8*1 .79,130-4*1 .79,3*2.82,100-3*2.82 
IF  d*4  then  7720 
GO  TO  7740 

VIEWPORT  8*1 .79,150-4*1 .79,3*2.82,100-3*2.82 

PRINT  90,17:1.2,1 .8 

CHARSIZE  1 

DELETE  K 

DIM  K(C0) 

CALL  "TIME",T$ 

PRINT  9D: "L";R$ 

PRINT  9D:Ti; 

REM 

WINDOW  LGT(N6) ,LGT(N7) ,LGT(N1 ) ,LGT(N2) 

IF  P>10  THEN  7980 
H1*INT(LGT(N1 >  > 

FOR  I*10^H1  TO  10*(H1  +  1)  STEP  10''H1*P 
MOVE  9D;LCT(N6) ,LCT(I) 


7860  DRAW  SD 5 LGT ( N7 > ,LGT ( I ) 

7870  NEXT  I 
7880 

7890  IF  Hl<aLGT<N2)  THEN  7840 
7900  HlsLGT(N6) 

7910  FOR  I=*tO''Ht  TO  10'^<Ht  +  l)  STEP  10^Hl*P 
7920  HOVE  eD;LGT<I) ,LGT(N1 ) 

7930  DRAW  SD : LGT < I ) , LGT ( N2 ) 

7940  NEXT  I 
7950  H1=:H1+1 

7960  IF  H1<=LGT<N7)  THEN  7910 

7970  GO  TO  8020 

7980  AXIS  @0; 1 ,1 ,0,-2 

7990  MOVE  eD;LGT<N6) ,LGT(N2) 

8000  DRAW  @D:LGT<N7) ,LGT(N2) 

8010  DRAW  eD:LGT<N7) ,LGT(N1 ) 

8020  PRINT  eD>17;2,3 
8030  CHARSIZE  4 

8040  GO  TO  R  OF  8050,8070,8090 
8050  A$=»MAXIMAX  SHOCK  SPECTRUM" 

8060  GO  TO  8100 

8070  A$= "RESIDUAL  SHOCK  SPECTRUM" 

8080  GO  TO  8100 

8090  Ai="MAXIMAX  and  RESIDUAL  SHOCK  SPECTRUM" 

8100  GOSUB  8170 

8110  IF  W1>0  THEN  8240 

8120  A$=W$ 

8130  PRINT  SD, 17; 1 .2,1 .8 
8140  CHARSIZE  2 
8150  GOSUB  8170 
8160  GO  TO  8240 

8170  MOVE  eD;LGT(10'^(0.5»LGT<N6itN7))),LGT(N2) 
8180  PRINT  eD: "KK"; 

8190  FOR  1=1  TO  LEN(At)/2 
8200  PRINT  8D: "H"; 

8210  NEXT  I 
8220  PRINT  eD;Ai; 

8230  RETURN 

8240  GO  TO  R  OF  8250,8270,8290 
8250  V$="P scud 0 velocity  n/sec " 

8260  GO  TO  8300 

8270  Vf=”Rel.  Velocity  «/sec" 

8280  GO  TO  8300 

8290  V$= "Vel oc i t  y  n/sec  " 

8300  MOVE  8D;LGT(N6>  >LGT(10*<0.5i>LGT(N2»Nl  )  )  ) 

8310  PRINT  eD,17; 1 .5,2.5 

8320  CHARSIZE  2 

8330  PRINT  SD; "HHHHHbb"; 

8340  FOR  1*1  TO  LEN(VS)/2 
8350  PRINT  SD; "K"; 

8360  NEXT  I 

8370  FOR  1=1  TO  LEN<V$) 

8380  X$=SEG(V$,I,1 ) 

8390  PRINT  eD;Xi; "HJ" j 
8400  NEXT  I 

8410  MOVE  eD;LGT( 10*<0 .5»LCT(N6*N7) ) > ,LGT(N1 ) 
8420  H$*"Ndtural  Frequency  Hz" 

8430  FOR  1*1  TO  LEN<H$)/2 
8440  PRINT  eD; "H”  j 
8450  NEXT  I 
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9460  PRINT  eO; 

8470  H1«0 

8480  FOR  I=L.GT<N1)  TO  LGT<N2>  STEP  lO^Hl 
8490  MOVE  eD:LGT<N6),I 

8500  PRINT  SD;  ■•HHHHH”;"10^";INT(H1+LGT(N1)) 

8510 

8520  NEXT  I 
8530  Hl^O 

8540  FOR  I=LGT(N6)  TO  LGT<N7)  STEP  lO^Hl 
8550  MOVE  eD! I,LGT(N1 ) 

8560  PRINT  SD: "HHJ” j "10^-;INT(H1+LGT(N6) ) 

8570  Hl=sHl  +  l 
8580  NEXT  I 

8590  MOVE  eD;LGT(N6) ,LGT<N1 > 

8600  IF  D*32  THEN  8620 
8610  PRINT  eD>17;l,2 
8620  CHARSIZE  2 
8630  A$»"(Zero  Daiiping)" 

8640  PRINT  eD;«JJJ”jA$ 

8650  GOSUB  8770 
8660  GOSUB  8820 
8670  IF  R=3  THEN  8690 
8680  RETURN 
8690  R=R>1 

8700  IF  D04  THEN  8740 

8710  PRINT  "LYOU  MAY  CHANGE  THE  PEN  COLOR  G"; 

8720  PRINT  "  Press  RETURN  to  continue” 

8730  INPUT  Ai 
8740  GOSUB  8770 
8750  GOSUB  8820 
8760  RETURN 

8770  GO  TO  R  OF  8780,8800,8780,8800 
8780  K=V5 
8790  GO  TO  8810 
8800  K=Y5 
8810  RETURN 
8820  1=1 

8830  FOR  I=I  TO  CO 

8840  IF  F5(I>>0.01  THEN  8870 

8850  NEXT  I 

8860  RETURN 

8870  MOVE  eD;LGT(F5(I) ) ,LCT<ABS<K<I) )  ) 

8880  FOR  I=I  TO  CO 
8890  IF  F5(I)=0.0l  THEN  8830 
8900  IF  R<*3  THEN  8920 
8910  DASH  170 

8920  DRAW  9D ; LGT < F5( I > ) ,LGT( ABS(K < I ) )  ) 

8930  NEXT  I 

8940  IF  W1 >0  THEN  8960 
8950  DASH  0 
8960  HOME 
8970  RETURN 

8980  REM  »****»*#  NUMERICAL  OUTPUT  ee**#e**ee*e»*eeee*e**eeeee*e* 

8990  WINDOW  0,130,0,100 

9000  VIEWPORT  0,130>0,100 

9010  I1»0 

9020  I2«l 

9030  N»1 

9040  IF  D*32  THEN  9070 
9050  IF  D«40  THEN  9130 
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9060  PRINT  @D>17:1 .5,2.5 

9070  CHARSIZE  3 

9080  PRINT 

9090  HOVE  80:65,99 

9100  FOR  I»1  TO  LEN(W^)/2 

9110  PRINT  90: 

9120  NEXT  I 
9130  PRINT  90:U$ 

9140  CHARSIZE  1 
9150  CALL  mHE”,T$ 

9160  PRINT  90:  jR»,  “JI” 

9170  IF  D-40  THEN  9190 

9180  PRINT  90:  USING  9190  ;  "KKK ;  "PAGE"  ;  12  j  " J J “ 

9190  IHAGE  3A,122X,4A,20,2A 
9200  CHARSIZE  2 

9210  PRINT  90:"  FREQ  SHOCK  SPECTRUM  FOURIER  TRANS’* 

9220  PRINT  90:'*  (HZ)  (MM)  (M/S)  (G)  (MS)  (M/S)  (OEC)^" 

9230  IMAGE  40 . 1 0, X , 50 . 20, X , 40 . 2D,X,5D . 10,X, 40 . 20,X, 40 . 20,2X, 40 
9240  IMAGE  65T , 40 . 1 0 ,X ,50 . 20 , X , 40 . 20 , X ,50 . ID , X , 40 . 20 ,X , 40 . 20 ,2X , 40 
9250  FOR  I-l+Il  TO  CO 
9260  IF  I>50»N  THEN  9300 

9270  PRINT  90;  USING  9230 : F5( I ) ,X5( I ) , V5( I > ,G5( I ) , T5 ( I ) , Y5( I ) , Z5( I ) 
9280  NEXT  I 
9290  RETURN 

9300  IF  O«40  THEN  9350 

9310  CHARSIZE  1 

9320  PRINT  8D:”^JJJJ1‘' 

9330  CHARSIZE  2 
9340  PRINT  90: "J" 

9350  FOR  I*l+50*N  TO  CO 
9360  IF  I>(N-H)«50  THEN  9430 
9370  IF  0*40  THEN  9400 

9380  PRINT  90:  USING  9240 : F5( I ) , X5( I ) , V5 ( I ) ,G5( I ) , T5( I ) , Y5( I ) , Z5( I ) 
9390  CO  TO  9410 

9400  PRINT  90;  USING  9230 ; F5 ( I ) , X5( I ) , V5 ( I ) ,G5( I ) ,T5( I ) , Y5( I ) , Z5( I ) 

9410  NEXT  I 

9420  RETURN 

9430  IF  0*40  THEN  9450 

9440  CALL  "WAIT", 4 

9450  11*100 

9460  I2*I2^l 

9470  N-N^2 

9480  IF  0*40  THEN  9500 
9490  GO  TO  9040 
9500  GO  TO  9250 

9510  REM  ANALYSIS 

9520  R*3 

9530  GOSUB  7670 

9540  R*1 

9550  MOVE  90:3.11,2 

9560  PRINT  90:  USING  "2A, 5X , 4A,3X ,7A" : "NR " , " CHZ 1 " , " CM/SECl " 

9570  POINTER  X,Y,Q8 

9580  IF  at*"a"  THEN  9750 

9590  MOVE  9D:X-LCT(1 . 022) ,Y+LCT(1 , 022) 

9600  CHARSIZE  2 

9610  PRINT  90: "J"";R 

9620  GOSUB  9650 

9630  R«R4>1 

9640  GO  TO  9570 

9650  VIEWPORT  106,130,0,100 
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9660  UINDOU  0,32>0,t00 
9670  MOVE  80:0,94 
9680  FOR  I»1  TO  R 
9690  PRINT 
9700  NEXT  I 

9710  PRINT  90:  USING  9720  :  " %R ,  ”  :  ”  ,  1 ,  1  O^Y 
9720  IMAGE  2A , 20 , A , 2X , 40 . 1 0 , 3X , 30 . 0 

9730  VIEWPORT  5^1  . 79 , 11  0-4)H  .  79 , 1 0+3»2 . 82 , 1  00-3^t2 . 82 
9740  WINDOW  LCT(N6> ,LGT(N7) ,LGT<N1 ) ,LGT(N2) 

9750  DASH  0 
9760  HOME 
9770  RETURN 

9780  REM  ik*^*^*it*****it***^**ii*  MULTI  DISPLAY 
9790  CHARSIZE  3 
9800  0>32 

9810  PRINT  -LTENTER  NUMBER  OF  DESIRED  SPECTRAG: 

9820  INPUT  wT 
9830  DIM  W2(W1) 

9840  PRINT  ".ENTER  FILE  NUMBERS  ON  DEVICE  NR . 3G :  ") 

9850  INPUT  U2 

9860  PRINT  “.ENTER;  1  FOR  MAXIMAX-SPECTRA.2  FOR  RESIDUAL-SPECTRAg :  " ; 
9870  INPUT  R 
9880  D3*W2(1) 

9890  COSUB  7570 

9900  GOSUB  7690 

9910  GOSUB  10060 

9920  IF  Wl*l  THEN  10050 

9930  FOR  10-2  TO  W1 

9940  D3-W2(I0> 

9950  COSUB  7570 
9960  DASH  11*10 
9970  GOSUB  8770 
9980  GOSUB  8820 
9990  GOSUB  10060 
10000  NEXT  10 
10010  DASH  0 
10020  Wl-0 

10030  DELETE  W2 , X , Y , 1 0 , D3 , Q$ 

10040  W6-"TEST" 

10050  RETURN 
10060  POINTER  X,Y,Q6 
10070  MOVE  9D:X,Y 
10080  DRAW  9D:X>0.25,Y 
10090  PRINT  BD; "  ">W$ 

10100  RETURN 
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1  REM  /////////////////  eSPRGM/PART  III  \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2  INIT 

3  GO  TO  260 

4  IF  0»3  THEN  94 

5  GO  TO  7450 

8  IF  0=*3  THEN  94 

9  GO  TO  7530 

12  IF  0=3  THEN  94 

13  GO  TO  8980 
16  GO  TO  190 
20  GO  TO  100 

24  IF  0=3  THEN  94 

25  R  =  1 

26  GO  TO  7690 

28  IF  0=3  THEN  94 

29  R*2 

30  GO  TO  7690 

32  IF  0=3  THEN  94 

33  Ra3 

34  GO  TO  7690 

36  IF  0=3  THEN  94 

37  GO  TO  9510 

40  IF  0=3  THEN  94 

41  CO  TO  9780 
45  RETURN 

48  IF  0=2  THEN  90 

49  GO  TO  510 

50  DELETE  A 

52  IF  0=2  THEN  90 

53  GO  TO  1780 

54  DELETE  A 

64  IF  0=3  THEN  94 

65  Rl*l 

66  GO  TO  6010 

68  IF  0=3  THEN  94 

69  Rl=4 

70  GO  TO  6010 

72  IF  0=3  THEN  94 

73  Rl*2 

74  GO  TO  6010 

76  IF  0=3  THEN  94 

77  Rl=3 

78  GO  TO  6010 

80  IF  0=3  THEN  94 

81  Rl=5 

82  GO  TO  6010 

90  PRINT  "You  got  to  PRESS  KEY  5  firstlG" 

92  RETURN 

94  PRINT  ’’You  got  to  PRESS  KEY  4  firstlG" 

96  RETURN 

100  load  ROUTINE 

110  IF  0=3  THEN  170 

120  DELETE  F5 , X5 ,U5 , C5 ,T5 , Y5 , 25 , C>C0 , K , A 
130  C0=90 
140  0=3 

150  APPEND  ’'9SPRGM/PARTIII"  ;  1 

160  PRINT  "JCSPRGM/PART  111  IS  LOADED!" 

170  DELETE  5140,20000 
180  RETURN 
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190  IF  0»2  THEN  240 
200  0»2 

210  DELETE  270>20000 

220  APPEND  "eSPRGM/PARTII‘*;260 

230  RENUMBER  6000 >10,260 

240  PRINT  -JGSPRGM/PART  II  IS  LOADED!" 

250  RETURN 

260  DEFAULT  CONDITION 

270  Wt="TEST" 

280  0-3 
290  Ul-0 
300  D-32 
310  Fl-1 
320  F2-1000 
330  F8-1000 
340  Nl-0.01 
350  N2-100 
360  N6-1 
370  N7-1000 
380  CO-90 
390  NO-4096 
400  Kl-2 
410  K2-2 
420  K3-10 
430  P-10 

440  R$-"NRL  •  CODE  5842" 

450  DIM  M9(7>2) >A2(7>2) >C(C0) 

460  M9-0 
470  A2-0 

480  DIM  DO (8) >Ht(20) >V0(22) >G$(7) ,M$(7>  >N$(8) >X$(22) >P*(35) >L$<4)>Qi<l ) 
490  PRINT  ”  GSPRCM/PART  III  IS  LOADED  (DEFAULT  CONDITIONS)!" 

500  RETURN 

510  REM  READ 

520  CHARSIZE  3 

530  PRINT  "LJENTER  FILE  NUMBER  (ON  DEVICE  3)  TO  READ  RECORD  FROMGj  "j 
540  INPUT  dT 
550  DIM  A(N0) 

560  U3-130 

570  IF  D04  THEN  600 
580  U3-150 

590  PRINT  8D>17j 1 .2>1 .8 
600  CHARSIZE  2 
610  CALL  "TIME">T$ 

620  PRINT  9Dt 
630  PRINT  9D:T$ 

640  FIND  B3!D1 

650  READ  93:U0>R2>H>T1 >M 

660  FIND  03:D1 

670  READ  93:W0>R2,H>T1 >M>M9(M>2) >A2(M>2) ,A 
680  M9<M,1)«M 
690  A2(M,1)«M 
700  S-M9<M,2> 

710  WINDOW  0>W3>0>100 
720  VIEWPORT  0>W3> 0,100 
730  IF  WDO  THEN  810 
740  PRINT  9D,17;2>3 
750  CHARSIZE  4 
760  MOVE  9D:W3/2>98 
770  FOR  I-l  TO  LEN<Wt)/2 
780  PRINT  9D:"H"; 
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790  NEXT  I 
800  PRINT  80:U$ 

810  PRINT  @0,17:1 .2,1 .8 
820  CHARSIZE  2 

830  GO  TO  R2  OF  840,970,1110 

840  GOSUB  1280 

850  S1»S 

860  GOSUB  1460 

870  PRINT  @0:”*" 

880  CALL  ”DIF3”,A,A 

890  S1=S1/H 

900  GOSUB  1340 

910  GOSUB  1460 

920  CALL  ”DIF3'*,A,A 

930  S1=S1»10^3/(H*9.8044) 

940  GOSUB  1400 
950  GOSUB  1460 
960  GO  TO  1230 
970  GOSUB  1340 
980  St*S 
990  GOSUB  1460 
1000  PRINT  8D: 

1010  CALL  -INT”,A,A 

1020  S1»S1)(H 

1030  GOSUB  1280 

1040  GOSUB  1460 

1050  CALL  ••DIF3“,A,A 

1060  CALL  **DIF3’*,A,A 

1070  Sl3Sl»10^3/(H^2^*9.8044) 

1080  GOSUB  1400 
1090  GOSUB  1460 
1100  GO  TO  1230 
1110  GOSUB  1400 
1120  S1»S 
1130  GOSUB  1460 
1140  PRINT  @D: 

1150  CALL  ”INT”,A,A 

1160  S1*SHH  . 0E“3*H*9.8044 

1170  GOSUB  1340 

1180  GOSUB  1460 

1190  CALL  '•INT",A,A 

1200  S1»S1»H 

1210  GOSUB  1280 

1220  GOSUB  1460 

1230  UIEWPORT  0,U3, 0,100 

1240  WINDOW  0,W3, 0,100 

1250  MOUE  90:0,0 

1260  PRINT  9D;*'»  orig.  tine  record"; 
1270  RETURN 

1200  VIEWPORT  5,W3-30,4,32 
1290  V^«"DISPLACEMENT" 

1300  G6>‘'D  max  :" 

1310  M$-"D  win  : " 

1320  Nf*"  [««!•* 

1330  RETURN 

1340  VIEWPORT  5 , W3-30 , 36 , 64 
1350  V$*"VELOCITY" 

1360  Gf>"V  wax 
1370  M$«"V  win 
1380  Cw/eecJ" 
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1390  RETURN 

1400  VIEWPORT  5,W3-30,68>96 
1410  V$»"ACCELERATION" 

1420  G««“A  nax  ; " 

1430  H$«”A  Min 
1440  Nfa”  Cq]” 

1450  RETURN 

1460  CALL  ”MAX%A,V1>I4 
1470  CALL  «MIN%A,V2>I3 
1480  WINDOW  0>N0>V2,V1 
1490  IF  D»32  THEN  1550 
1500  HOVE  8D:0,0 
1510  FOR  lal  TO  NO 
1520  DRAW  8D:I,A(I) 

1530  NEXT  I 

1540  GO  TO  1560 

1550  CALL  "DISP^A 

1560  AXIS  8Dt400,0 

1570  FOR  I»0  TO  NO  STEP  2000 

1580  HOVE  eD:I,0 

1590  PRINT  eD:"HJ",I*H; 

1600  NEXT  I 

1610  PRINT  SD: -HHHHHHHHHHHCMS]" 

1620  MOVE  8D:4200>vT 

1630  PRINT  8D:  USING  1670 : ,G$>V1«S1 ,Ni; 

1640  MOVE  eD; 4200, VI 
1650  PRINT  eD! j 

1660  PRINT  SD;  USING  1 670 : ” J ” , M$ , V2»S1 , N$ ; 

1670  IMAGE  A>7A,6D. 1D>11A 
1680  MOVE  8Ds0>0.5«(V2'(-Vl) 

1690  PRINT  8D:**HH»; 

1700  FOR  I»1  TO  LEN<V$)/2 
1710  PRINT  8D: "K"; 

1720  NEXT  I 

1730  FOR  I«1  TO  LEN(V$) 

1740  X$»SEG(Vi,I,l  ) 

1750  PRINT  8DjX$j "HJ"; 

1760  NEXT  I 
1770  RETURN 

1780  REM  PROCESS 

1790  CHARSIZE  3 
1800  CALL  mME**>Bi 
1810  PRINT  8Di  “i’* 

1820  PRINT  832,26:2 
1830  IF  F1<1  THEN  1890 
1840  IF  FDIOOO  THEN  1890 
1850  IF  F2<1  THEN  1900 
1860  IF  F2>1000  THEN  1900 
1870  IF  F2<*F1  THEN  1900 
1880  GO  TO  1930 
1890  IF  C0«1  THEN  1930 

1900  PRI  -ANALYSIS  FREQUENCY  IS  OUT  OF  RANGE  ENTER  NEW  VALUES  OF  F1,F2G“ 
1910  INPUT  F1,F2 
1920  GO  TO  1800 

1930  DELETE  F5,X5,V5,G5,T5, Y5,Z5,C,K ,A 
1940  DIM  A(N0) 

1950  C3«INT< (MEM0RY/8-18)/a> 

1960  DIM  F5<C0> ,X5<C0) ,V5(C0) ,G5<C0) ,T5<C0) ,Y5(C0> ,Z5(C0) ,C(C0> 

1970  FIND  83:D1 

1980  READ  83:W6,R2>H,T1 ,M,S 
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1990  FIND  @3:01 

2000  E1=INT(1000/(K1»F8«H)) 

2010  E2*INT( 1000/(K2»F8»H) ) 

2020  E3=INT(1000/(K3*F8»H) ) 

2030  F5( 1 )=F1 

2040  IF  C0=1  THEN  2080 

2050  FOR  1=2  TO  CO 

2060  F5(I)=(F2/F1 )‘(1/<C0-1 ) )xF5<I-l > 

2070  NEXT  I 
2080  X5=0.01 
2090  V5=0.01 
2100  05=0.01 
2110  T5=0.01 
2120  Y5=0.01 
2130  Z5=0.01 
2140  11=1 
2150  12=0 

2160  IF  Fl<10  THEN  2200 
2170  12=1 

2180  IF  Fl<100  THEN  2200 

2190  12=2 

2200  GO  TO  2370 

2210  E=E1 

2220  GOSUB  2460 

2230  IF  I1>C0  THEN  2310 

2240  GO  TO  2370 

'3p^,n  FsE2 

2260  GOSUB  2460 

2270  IF  I1>C0  THEN  2310 

2280  GO  TO  2370 

2290  E=E3 

2300  GOSUB  2460 

2310  DELETE  A 

2320  CALL  "TIME",T$ 

2330  GOSUB  3320 
2340  PRINT  @32,26:0 

2350  PRINT  eD:"_MAX.  POSSIBLE  ANALYSIS  FREQUENCIES:  ;C3 
2360  RETURN 
2370  F3=F5<I1) 

2380  C1=I1 

2390  12=12+1 

2400  FOR  11=1  TO  CO 

2410  IF  F5(ll>>10"12  THEN  2430 

2420  NEXT  II 

2430  F4=F5< Il-l  ) 

2440  CO  TO  12  OF  2210,2250,2290 
2450  RETURN 

2460  READ  @3 : W* , R2 , H , T1 , M , S , A2(M , 2 > , A 

2470  FIND  e3:Dl 

2480  IF  E<=1  THEN  2550 

2490  J=1 

2500  FOR  1=1  TO  NO  STEP  E 
2510  A(J)=A<I> 

2520  J=J+1 
2530  NEXT  I 
2540  GO  TO  2560 
2550  E=1 
2560  N5»N0/E 
2570  H2»H»E 

2580  GO  TO  R2  OF  2590,2630,2660 
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2590  CALL  “DIF3“,A>A 

2600  CALL  ‘•DIF3”,A,A 

2610  S2*S 

2620  CO  TO  2670 

2630  CALL  "DIF3",A,A 

2640  S2*S»H2 

2650  GO  TO  2670 

2660  S2=S»10^-3^tH2^2»9.8044 

2670  Ca*!!-! 

2680  PI-0 . 0062831853»F3)(H2 
2690  FOR  I=C1  TO  C2 
2700  C(I)=2»C0S<P1 ) 

2710  IF  C1=C2  THEN  2730 
2720  P1=P1)KF4/F3>^(1/(C2“C1)) 

2730  NEXT  I 

2740  IF  Pl<0.902  THEN  2790 
2750  IF  143.566/H2>F2  THEN  2790 

2760  PRINT  ’'LGWARNING;  FREQUENCIES  ARE  TOO  HIGH  FOR  YOUR  TIME  INCREMENT" 
2770  PRINT  '*  ERRORS  ARE  OUTSIDE  THE  RANGE  -7  TO  +3  PERCENT  FOR" 

2780  PRINT  "  FREQUENCIES  ABOVE  " ) 143 , 566/H2 j "  HERTZ" 

2790  FOR  I=C1  TO  C2 
2800  F=F5(I) 

2810  IF  F<-1  THEN  2840 
2820  Nsl000/(H2*5*LGT(F) ) 

2830  IF  N<N5  THEN  2850 

2840  N=N5 

2850  X1=0 

2860  X2=0 

28’0  R3*0 

2880  R4=0 

2890  IF  F=0  THEN  2920 
2900  N4=N-l-500/<F)*H2) 

2910  IF  N4>3  THEN  2930 
2920  N4*3 

2930  FOR  J*1  TO  N-1 
2940  X3*C(I)*X2-Xl-»'A(  J) 

2950  IF  J<N4  THEN  3000 

2960  IF  ABS(X3)<ABS<R4)  THEN  3030 

2970  R4aX3 

2980  T4*J 

2990  GO  TO  3030 

3000  IF  ABS(X3)<ABS<R3)  THEN  3030 

3010  R3=X3 

3020  T3*J 

3030  X1=X2 

3040  X2*X3 

3050  NEXT  J 

3060  IF  ABS<R3) >ABS(R4)  THEN  3090 

3070  R3=R4 

3080  T3==T4 

3090  X*S2«tR3 

3100  V^O  .  006283184(F«X 

3110  G=6.288318»Fi^V/9.8044 

3120  T*Tl+(T3-0 .5)«H2 

3130  Y=:0.00628318»F»S2»R4 

3140  Z=1  .  OE-3^KT1  +  (T4-0 .5)  xH2)«F+0 .25 

3150  Z*Z-INT(Z) 

3160  IF  Y>0  THEN  3180 
3170  2*Z-0.S 

3180  IF  Z<*0 .5  THEN  3200 
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3190  ZsZ-1 
3200  Z=360*Z 
3210  F5(I)^F 
3220  X5<I)aX 
3230  y5(I)=»V 
3240  G5<I)aG 
3250  T5(I)«T 
3260  Y5(I)=Y 
3270  Z5(I)=Z 

3280  PRINT  eO:  USING  3290 : I ,F>U,N; 

3290  IMAGE  20T , 3D, 4X , 4D . 1 D , 4X , 5D . 2D, 3X ,6D 
3300  NEXT  I 
3310  RETURN 

3320  REM  **!****  SUBROUTINE  TIME 
3330  DIM  B(2) 

3340 

3350  FOR  lal  TO  2 
3360  Xf=»SEG<CS,ll  ,8) 

3370  Y$*SEG<X$,4,2) 

3380  Zi*SEG<X$,7,2) 

3390  B<I )«UAL(X$)+(VAL< Y$)»60+VAL(2$) >/3600 
3400  C$*Ti 
3410  NEXT  I 
3420  U*B(2)-B(1) 

3430  B3a(U-INT<U) )*60 
3440  B4®<B3'-INT(B3)  )  »60 

3450  PRINT  ©D;  “^CALCULATION  TIME  (hM  Min:  sec); 

3460  PRINT  8D:  USING  3470  ;  INT<U)  ,  ”  :  ”  ,  INT<B3)  ,  "  :  •MNT(B4) 

3470  IMAGE  2D, A, 2D, A, 2D 
3480  RETURN 

5000  REM - - - 

5010  REM - - 

5020  REM  THIS  PROCEDURE  MAY  BE  USED  TO  CALCULATE  SHOCK  SPECTRA  OF 

5030  REM  TIME  RECORDS  STORED  ON  EXT.  TAPE  (FILE  1  TO  7)  AND  TO  STORE 

5040  REM  THE  ACCORDING  DATA  ON  INT.  TAPE  (FILE  TO  FILE  »+7). 

5050  FOR  Dl«l  TO  7 
5060  GOSUB  52 
5070  N=Dl+0 
5080  FIND  N 

5090  MARK  1 ,7000 

5100  FIND  N 

5110  WRITE  C0,Wi,F5,X5,U5,G5,T5,Y5,Z5 
5120  NEXT  D1 
5130  RETURN 
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